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ABSTRACT 


Most treat ments of the problem of temperature distribution in a solid body during 
heating or cooling assume that the change in surface temperature is either instan- 
taneous or linear with time. Experience shows that in many cases of practical impor- 
tance the increase in surface temperature may be represented by a simple exponential 
function of the time. Equations for the distribution of temperature within a body when 
the surface is heated in accordance with this exponential relation have been derived for 
the following shapes: (1) rectangular parallelopiped (brick); (1a) rectangular rod; 
(1b) slab; (2) cylinder: (2a) cylindrical rod; (3) sphere. These equations also represent 
the general case of diffusion under the given boundary conditions and apply, therefore, 
equally well to problems such as the diffusion of one substance into another, as for 
instance, one metal into another. 


INTRODUCTION 


HE problem of the distribution of temperature within a solid body of 

simple geometrical shape has been discussed many times but the treat- 
ment has almost invariably been based on an assumption either of a linear 
rate of heating of the surface of the body or of an instantaneous change in the 
skin temperature from one constant value to another. The equations so de- 
rived, while quite appropriate for certain purposes, cannot, without modifica- 
tion, be applied to many processes of practical importance, such as the heating 
of cold metal masses by immersion in a hot furnace, when as is well known, 
the change of surface temperature cannot be regarded as following either of 
these assumptions. 

Experimental work by the author on small steel cylinders introduced into 
a hot furnace, or heated from room temperature with the furnace, indicates 
that the surface temperature may be represented by an exponential function 
of the form 6 = 6’(1 —e-**) where @ is the surface temperature at time ¢ and 0’ 
is the final equilibrium temperature of the surface.* The problem then, is to 
find equations expressing the temperature distribution within a solid when 
the surface is heated in accordance with this relation. The equations which 
have been derived as solutions of this particular problem are, of course, solu- 
tions of the general problem of diffusion in solid bodies under the aforemen- 


* If the surface is cooled the exponential function is = 0’e~**, 
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tioned boundary conditions; they are not limited to a description of tempera- 
ture distribution but may be applied to other problems of such as the diffusion 
of one substance into another, as for instance, one metal into another, pro- 
vided always that the surface concentration varies in the manner represented 
by the exponential function. In these cases it is necessary only to let the 
symbols for temperature represent the corresponding concentration of the 
diffusing substance. 

Since the application in which the author is chiefly interested is to the 
heat treatment of metals and since the important item, in this case, is the 
temperature distribution near the center, the solution has been given in terms 
of sine or cosine series since they are more convergent in this region than 
the probability integral, which is an alternative form of solution. The equa- 
tions are, however, rigorous and perfectly general; they apply to any portion 
of the solid body or may be used in other ways. 

For instance, anticipating a little, one of the Eqs. (41) indicates a rela- 
tively simple method for the determination of the heat diffusivity of solids. 
This equation gives an expression for the change of temperature at the center 
of a cylinder, initially at uniform temperature, when placed in a furnace at 
an elevated temperature. The necessary measurements can be made by means 
of a thermocouple attached to, or just within, the surface and another placed 
in a small axial hole. The latter will not appreciably disturb the equilibrium 
at the center since the gradient across it is zero. 

The general method is that of Williamson and Adams! who considered 
the cases of linear and instantaneous change of surface temperature. Since 
their paper contains several typographical errors their method has been re- 
produced in some detail in the present paper in order that the procedure may 
be made clear. 

Solids of the following shapes are considered: (1) Brick (rectangular paral- 
lelopiped); (1a) Rod of rectangular section, or brick with two opposite faces 
insulated (infinite rod); (1b) Infinite slab or thin plate; (2) Cylinder; (2a) 
Cylindrical rod or cylinder with insulated ends (infinite cylinder) ; (3) Sphere. 
In all cases the initial temperature is taken as uniform.* 


NOMENCLATURE 

6 =difference from the initial temperature? at the point (x, y, z) at time 
t seconds from start. We shall in all cases consider the initial tem- 
perature to be 0 so that # represents the temperature at time ¢ seconds 
later. 

«x =thermal diffusivity in cm* per second 

6’ = final temperature at the surface of the metal, that is, the temperature 
of furnace 

t= time in seconds 


1 Williamson and Adams, Phys. Rev. 14, 99 (1919). 

? For problems of diffusion of one substance into another substitute the appropriate con- 
centration for temperature throughout. 

* The equations as given have been derived for the case of heating of the surface. To apply 
them to problems of cooling 6=0’(1—e~**) is replaced by 0e~*' in the suggesged form of solu- 
tion. This requires a similar substitution in the other equations and changes the sign of the 
summation term. 
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TEMPERATURE DISTRIBUTION IN SOLIDS 


Williamson and Adams have shown that a solution of the differential 
equations for heat conduction or diffusion may be made up of terms which are 
the product of an exponential function of the time and a trigonometrical func- 
tion of the coordinates, except in case (2) where Bessel’s functions are re- 
quired. The general form of the solution is 

marx | Ty pr: 


- sin —— sin ——(1 — e**), (1) 
a b c 


6 = 6'(1 — 8) + Sof(m, n, p) sin 





This satisfies the boundary conditions that @ must vanish for f=0 and must 
approach @ as ¢ increases. The exponent 8 must be determined experimentally 
but we may determine f(m, 1, p) and a@ analytically, using the general equa- 
tions for heat conduction. Thus for the brick we have 


06 0°76 0°60 0°60 

—=k a marae: a nen 

ot Ox dy? 02? 
By differentiating (1) we may obtain the values of the various differential 


coefficients in (2) and thus find values for a and f (m,n, p) which satisfy the 
above conditions. A number of specific cases are now considered. 





Case 1. Brick shape (rectangular parallelopiped). Origin of coordinates 
at thecenter of the block. Boundary conditions, @ = 0, wheret = 0,6 =6'(1—e **) 
over the surface x = +a; y= +);s= +c. 

We have then: 


mr x nny prz 
0 = 0'(1 — e#) + > f(m, n, p) cos —— cos —— cos ——(1 — e~**) (3) 
2a 2b 2c 
where m, n, p are integers. The equation of conduction is 
6 
— = kA’ (4) 
at 


Following our general method we differentiate (3) and put the coefficients 
in (4) 


, mT Xx nr y prz 
6’Be-8t + So ae-#tf(m, n, p) cos ZT cos = cos 
2a 


mr? nx? p*n? 
=—k ) f(m, n, (“ + — + — -) 
DIK p) 4a? 4° 4c? 


2c 


mr Xx nury prz ° 
cos cos cos ——(1 — e") (5) 
2a 2b 2¢ 





Each side of this identity is made up of t-vo parts, one containing e-*! 
and the other not. Equating these separately yields: 


MT Xx nny prz 
cos - COS 
2c 











dae-**f(m, n, p) cos 
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12a? 2-2 2-2 
— k Dof(m, n, p(— + nit + “~) 











45? 4c? 
NX unr y prz 
cos ——— cos —— cos —— e~*! (6) 
2a 2b 2c 
and 
mr? ne? pr? \ mmx nr y pre 
6'Be-8t = — k > f(m, n, pf) (“ + —— + —— } 008 ———— cos —— cos ——— (7) 
4a?* 45° 4c* 2a 26 2¢ 


From (6) we obtain at once 


mr? nar? px’ 
leant. + + ‘~) (8) 
4a? 4b? 4¢? 


To evaluate (7) we need the series 





WX 1 Sax 1 Srx . 
—S Oi - — = > — os ——— > - <> (9) 
4 2a 3 2a 5 2a 
which holds between x = +a and x= —a. 


Multiplying together 3 such series yields 


= D(a 1yntete 


cos[(2m—1)x )arx/ (2a }cos[(2n—1)xy/ 2b |cos|(2p—1) mz, 2c] 
7 i 


10) 
(2m — 1)(2n — 1)(2p — 1) 


where m, n, p are all positive integers. 
Eq. (7) will reduce to this if f(m, 7, p) be put equal to 





640’ Be (— 1) m+n+ptl 
—_  _————————————————— ——_——— (11) 
kr* { (2m —1)?? (2n—1)?x? (2p—1)?x*) 
en Ee —_—_—_——_ + ———_ ———— 
4a? 4b? 42S 
Substitution of these values in (3) gives: 
640’Be8t P.n_m=x 
= #(1 — e 4) — —— dD (= 1) mtnter 
kr® mn, p=1 
cos [(2m — I)rx, 2a] cos [(2n — I)ry, 2b] cos oe. - AD) )arz 2c| ( 
Eee ~ — (12) 
j ( (2m — 1)? (2n — 1)% (2p — 1)*x?) 
(2m — 1)(2n — 1)(2p — 1)4 —_——— +- ——— 
4 4a? 4b? 42 
(2m — 1)?x* (2n—1 (2 ae 
(1 — exp — Sasol ‘a> om te + = =t/) 
(4a? 48° 42S 


By substituting Qn, Qn, Qp for (2m—1)x/2, (2n—1)m/2, (2b—1)x/2 re- 
spectively we get 
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8@’ eB p,n,m=o 
6 = O(1 — eS) — Sepe™ pm (— 1) materi 


hk m,n, p=1 
cos Q,x/a cos Qny/b cos 03/4 ‘Cc 
P* 


Dns On, On(Om 2/a? a QO,” / 52 — 0,?, ) 2/ jay’ 





1 — e~ k1Qm? a2+Q,,2/ 62+ Qp?/<?} ¢') (13) 


Case 1 (a). Rectangular rod. This case may be deduced directly from(13) 
by assuming that c is large compared to a and b. 


40’ Bet m ,n=0 


A= (1 — e~ Bt) _ p (— 1) *ta+1 
k m n=l 
cos (Qma/a) cos (( en — ek 1Qn7/a7+Q," ‘/ b*) ) (14) 
Om; Qn(On?/a® + Qu? | 


Case 1 (b). Infinite slab. This case also may be deduced from (13) by 
assuming both } and c as large compared to a. 





20’ a amet cos (Onx«/a) 
= 6’(1 —e aad ae ~ 1) +1. wth ee (1 Eyes e~ *(Qm?/a2) t) (15) 
m=1 On?) a“ 


Case 2. Cylinder. Radius =a, length = 2b. Origin of coordinates at center. 
In this case the computation is simplified if we introduce Bessel’s functions 
in place of the cosine terms in the radial direction. Only the following proper- 
ties of Bessel’s functions are required in the derivation. 


-2 6 











x z* x 
J(x) =1-—-—+— a eee (16) 
2? 2? xX 4 2? XK 4? X 6? 
dJ o(x) 
= —J;(x) (17) 
dx 
dJ ;(x) 1 
— = Jo(x) — —Ji(x) (18) 
dx x 
m=eo 2 x 
p a am ® “) = 1, (19) 
m=1 R,J (Rn) a 
The suggested form of solution is: 
. x mmZ 
= 6(1—e*) + Yy(m, RI Re =) cos (1 — e-*#), (20) 
a é 
The equation of conduction is 
rele 0? 1 086 076 
— se (— ra 2 =): (21) 
ot Ox* Z 0% 02° 


By reasoning similar to that used in considering the brick and by using 
(9) and (19) we obtain 
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80’ Be n,m=0 (— 1)™+! 
= 6(1 — eA) — — 








wk m n=1 (2m a 1) 
Jo(R,(x/a)) cos ((2m — 1)xz/2b)) ( R,,? (2m — 1)?x? 
_— ——— {1 —exp— Rd ——— hy) (22) 
R,2 (2m — 1)?x? la 4h? f 
( + - )RuI(R,) 
a? 4b? 


Case 2 (a). Cylindrical rod. When ) becomes large compared to a we get 





Oa V1 ~ %) —- —— —exp—k—t 


VA! Rp—Bty2 n=~ . > 2 
20’ Be a? abn (1 R,, ), (23) 
k n=1 R,°J \(R») a“ 


Case 3. Sphere of radius a, origin at center. The suggested form of solu- 
tion is 


m= 


6=6'(1 — e*) + > f(m) 
m= 1 


1 mr x 
da ——{) + o*). (24) 


x a 
The equation of conduction is: 
0( x6) ; 9?( x6) 
Ot Ox? 





(25) 


Following our usual procedure we get 








_ mnx mr? marx 
x0’BeF* +- > ae-**f(m) sin —— = — k )f(m) —ma—{i—- oo") (26) 
a a°* a 
whence 
9 9 
m7? _ 
qa — §— (27) 
a’ 
and since 
x 2a 
>——(- 1)" = — (28) 


- O’Be-8* 2a a? 
*. f(m) = — “neewrwe (29) 
k mar m1 


Hence 


eC ic 8 = (— 1) «¢ mrs ae 
O= (1 =) — ———. Fs —_— — sin —01 — i mtetietot), (30) 
k m=1 m3 x a 
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In many cases we are interested in the temperature at the center of the 
body (x =0, v=0, s=0). In this case all the cosine terms and J) terms become 
unity and we obtain 


Brick 
S6’Be~8! Pn mae 
6 = 6'(1 — et) — > (= A) tntert 
k mien Pp 1 
1 GC, 0,* Q 
$$$ $$$ (1 — exp — k & +> — ) (31) 
) QO (+ “-.) \ a* b? ¢* 
-” ( " . eee t — 
. a? 2 c2 


Rectangular rod 
40’Be Bt nu m=x 


g=0'(1—e  Jaaaoee S } & (— 1)mtntt 
4 m .n=1 


1 Om" Qn*) ’ 
a (1 —exp— ki" 4+~—* ‘) (32) 
0,,* 0, \ a~ b? f 


On; On; ( ri + SS ) 
a* b? 


Infinite slab 


— (1 — o-6bt _ 20'Be** " an ake a —iillle Om" , 
6= 0'(1 — e**) > ( 1) 1 — exp k i}. (33) 








k m=1 0." a’ 
Cylinder 
80’Be Bt n,m=2 
oO — ee) ES (ym 
wk m n=1 
1 j &,? ma?) 
ee ep 2 + - 1). (34) 
(2 1) (“ 4 =) RAR.) l a pz f 
m — — fp ee | R,, 
a? 4b? 


Cylindrical rod 


20a70'Be8* m=* 1 R,? 
6= 6(1 — eft) — ——_—_— a (1 —exp — k - -{ ) (35) 
k a R,°J (CR, ) : 


a? 
Sphere 
2a76’Be-Bt m=" (— 1) mr? 
6= 60'(1 — eft) — ———_—. »z — (1 —exp—k- 1), (36) 
k m-=1 ni? a? 


A still further simplification may be introduced. As the time ¢ increases 
the exponential terms approach zero and the terms involving the product 
of e~*‘ and e-** in the final equations are small compared to the others. This 
simplification may be applied to any of the equations but the ones for the 
center of the bodies are given below: 
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Brick 
80’BeBt Pn mmee (~ 1) mtntptt 
6= 6'(1 — et) — = > _—__ an 0 e" 0 sane, (37) 
m np=1 m n Pp 
Om, Qn, Qo, (Sv + a pe es ~) 
a b? c 
Rectangular rod 
46'Be-** n,m =0 _ ]) *+et! 
= §’(1 — e Ft) — — — ————_———, 38 
( hk m,n ,=1 On? G,* ( 
Va, On; (= ae + =] 
a* b? 
Infinite slab 
2476’ Be Bt m=n0 ais 1 m+l1 
@ = O(1 — oH) — —— 2. (39) 
h m=1 0,° 
Cylinder 
80’ Bet n=2 m= * 1)™+! 
iaN-¢)<—— Fo... (40) 
wk bined aeeall R,,? mm 
(2m — 1) eS — + - —-) RJ (RK) 
Cylindrical rod 
20’Be8! n=” 1 
6=0'(1 — e%) — ——_——_ ——--—-—*, 41 
k X RFI (Rn) (41) 
Sphere 
20’ Be- Bt m=2 =— 1) +! 
6= 0'(1 — eft) — — ee . (42) 





k a=-1 mx? 


Eqs. (37) to (42) have the advantage that they make it possible to ex- 
press k explicitly in terms of the other variables. In all the relations previ- 
ously discussed this is not easily accomplished, if indeed it is possible at all, 
and if k is to be determined resort must be had to some method of approxima- 
tion. The foregoing equations and the ones to follow are thus especially 
adapted to the measurement of the diffusivity. 


It will also be observed that in the equations for the so-called “unidimen- 
sional solids,” i.e., solids in which heat flows in one dimension only as for in- 
stance in the slab (39), cylindrical rod (41) and sphere (42), the term under 
the summation sign does not involve any of the variables in the problem; 
hence, remembering that 0, =(2m—1)2/2 we may replace the summation 
term by a numerical value. This operation leads to the following equations 
which have been rearranged to show their similarity of form: 
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Infinite slab 


2a70’Be#t 








6= 6(1 — e Ft) — — (0.2308). (43) 
Cylindrical rod 
2070’ Be-## 
6= 0 (1 — e*) — ——_——-(0. 1247). (44) 
Sphere 
2a°6’Be-#" 
6 = 6'(1 — et) — —— (0.0291). (45) 


The equations for the cube and for the rod of square cross section, can, 
obviously, be made to assume the same form. 
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Fig. 1. Temperature distribution in solid bodies during heating or cooling. 


The distribution of temperature within the body described by these equa- 
tions is quite similar to that obtained with linear heating of the surface; such 
a distribution is well known and has been shown by many others, so that no 
discussion need be given here. It is interesting, however, to investigate the 
variation with time of the temperature at the center of bodies of simple 
geometrical shape under the heating conditions postulated. This variation 
has been calculated for the unidimensional solids and is shown in Fig. 1. The 
following values have been assumed in making the calculations 6’ = 1000°C. 
B=2X10-, k=0.010 cm?/sec.,a=7.5cm. The values of 8 and k are rounded 
values based on on actual experience of the author in dealing with steel cylin- 
ders. 

It will be seen at once, that for a given thickness, or radius, temperature 
equalization is attained most rapidly in a sphere and least rapidly in a slab, 
which is the order observed with the common boundary conditions. Other 
general relations of the curves of Fig. 1 are obvious and need no further dis- 
cussion. These remarks apply, mutatis mutandis, to other cases of diffusion 
under the same boundary conditions. 
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ABSTRACT 

An investigation has been made of some thermionic properties of a barium film 
adsorbed on tungsten for various values of 0, the fraction of the surface covered with 
barium. The experimental tube was so constructed that the measurements could be 
made on a small area of the filament in order that errors due to lead losses, non-uni- 
formities, and voltage drop along the filament might be minimized. An electrometer 
vacuum tube was used to measure the electron currents. The experimental results ob- 
tained by Eglin were confirmed, i.e. the thermionic activity of the composite surface 
was found to increase with @ until an optimum value was reached after which further 
deposit of barium caused the activity to decrease. Asa result of prolonged deposition, 
a steady value of the activity was reached which did not change appreciably or regu- 
larly with further deposition. The change in C.D.P. was found to be directly propor- 
tional to the change in 7%, the absolute temperature of the filament required to give 
an arbitrarily fixed value of electron emission. It is significant that the factor of propor- 
tionality for all values of the activity higher than the steady value reached after 
prolonged deposition was found to be different than that obtained for the values 
which were lower. The distribution of velocities of the emitted electrons was found to 
be Maxwellian for all values of @. The construction of the experimental tube was found 
to be unsuitable for studying the effect of field intensity on the activity. Thus the 
results obtained pertaining to this effect are of little significance. 


T 1S a well-known fact that the theories of thermionic emission as derived 
for clean metal surfaces do not fit experimental facts obtained for com- 
posite surfaces. The coefficient A of the thermionic emission equation 
I=AT*e—¢/*? which theoretically is a universal constant for pure metals 
varies widely for composite surfaces. The Schottky theory accounts satis- 
factorily for the effect of field strength on the emission from clean metal sur- 
faces. It fails, however, to account for the effects observed in the case of com- 
posite surfaces. Experimental evidence indicates that @ for clean metals is 
independent of temperature. The results obtained for composite surfaces, on 
the other hand, have led J. A. Becker,! Schottky,? and others to the belief 
that the work function of a composite surface is dependent upon the tem- 
perature. 
Theories have been advanced by Langmuir,’ Becker,‘ and others to ac- 
count for the thermionic behavior of composite surfaces. Though these theo- 


1 J. A. Becker, Phys. Rev. 34, 1323 (1929). 

2 Schottky, Rothe and Simon, vol. 13, of Wien Harms, “Handbuch der Experimental- 
physic.” 

3 Langmuir, Gen. Elec. Rev. 23, 504 (1920). 

* Becker and Mueller, Phys. Rev. 31, 431 (1928). 
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ries may be qualitatively correct they are apparently not complete in so far 
as they fail to account for the effect of field strength on the emission at low 
accelerating potentials.® 

In view of the above, an investigation which might throw more light on 
the thermionics of composite surfaces appeared to be in order. In this in- 
vestigation, therefore, the following thermionic phenomena were studied: I. 
The change in thermionic activity and contact potential difference (C.D.P.) 
caused by depositing barium on tungsten. II. The relation between change 
in activity and change in C.D.P. III. The effect of field strength on the activ- 
ity and on the C.D.P. IV. The distribution of velocities of the emitted elec- 
trons for various values of @, the fraction of the surface covered with barium. 


EXPERIMENTAL ARRANGEMENT 


A diagram of the apparatus is shown in Fig. 1. The experimental tube 
contains the anode A, the collector C, the tungsten ribbon filament W, and 
the barium oxide coated filaments O; and O,. The anode consists of two paral- 
lel plates electrically connected. The filaments are kept taut by means of 
tungsten springs. The oxide coated filaments are mounted perpendicularly 
with respect to the tungsten filament in such a manner that the “beams” of 




















=S==— + 
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Fig. 1. Diagram of apparatus. 


material evaporated from the barium oxide coating are “focused” on the 
point P of the tungsten ribbon directly opposite the aperture // in the anode. 
(In this experiment the filaments were glowed under such conditions that the 
material evaporated could be assumed to be barium.)! The plate C serves 
as a collector for the electron beam through the aperture. 

The anode plates were 7.5 cm long, 3.0 cm wide, and were spaced 1.2 cm 
apart. The average thickness and width of the tungsten filament W were 
0.0025 cm and 0.0497 cm respectively. It was 7.0 cm long. The aperture // 
was 0.1 cm in diameter. The perpendicular distance between the aperture J/ 
and the tungsten ribbon was 0.6 cm. 

Precautions were taken to insure a high leakage resistance between the 
collector and the other electrodes of the tube. The electron current to C, J;, 
was measured by means of the electrometer tube arrangement £.° The electro- 
meter tube and the experimental tube were enclosed in a brass tube closed 


5 Compton and Langmuir, Reviews of Modern Physics 2, 150 (1930). 
* H. Nelson, The Rev. of Scientific Instr. 1, 281 (1930). 
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at the ends by bakelite disks through which the leads were brought out. The 
brass tube served as an electrostatic shield and as a means of keeping the 
tubes in dry air to prevent electrical leakage over the glass. 

The experimental tubes used in the investigation were thoroughly ex- 
hausted before being sealed off. The tubes were baked during exhaust and the 
metal parts were degassed. After the tubes had been sealed off from the pump- 
ing system, the vacuum was further improved by the glowing of auxiliary 
thoriated tungsten filaments. The excellent reproducibility of experimental 
results indicated that good vacuum conditions were obtained. 


MetTHOD 


A relative measure of the activity of the tungsten surface for various 
degrees of covering with barium was obtained by determining the absolute 
temperature 7 of the filament necessary to cause an electron current of 
1.9X10-” amperes to flow to the collector C, the anode and the collector 
being fixed at a positive potential of six volts with respect to the point P on 
the tungsten ribbon. 

From the Richardson equation of thermionic emission 


I = AT%e~¢elkT 


it follows that the relation between 7°, and the commonly accepted measure of 
thermionic activity, ¢, may be expressed as follows: 


oo = T,/11600[In a + 2 In Ty — In (1.9 X 10-12) + In Ao] 


where ¢¢ and A, are the values of ¢ and A, expressed in volts and amp./cm? 
deg’. respectively, when the fraction of the surface covered is equal to @ and 
where a is the effective emitting area in cm?. The value 11,600 deg./volt has 
been substituted for the ratio e/k. 

Though the value of a may be obtained experimentally by measuring the 
electron current to the collector at a known temperature of the clean tungsten 
ribbon, it is apparent that the value of ¢s may be determined through meas- 
urements of 7 only when Ag is known. 

The temperature at the point P of the tungsten filament corresponding to 
a definite filament current was determined from calculations based upon a 
knowledge of the dimensions of the ribbon and data published by Langmuir 
and Jones.’ End loss corrections were calculated in accodrance with the best 
published data.* An experimental check on the calculated results was ob- 
tained through pyrometer measurements. In Fig. 2 the temperature of point 
P of the tungsten ribbon is plotted as a function of the filament current. The 
points marked by circles are the calculated values; those marked by triangles 
are the ones obtained through pyrometer measurements. The points marked 
by crosses were calculated from Maxwellian distribution curves. 

The contact potential difference, V., between the point P on the tungsten 
ribbon and the collector C was determined from Maxwellian distribution 


7H. A. Jones and I. Langmuir, Gen. Elec. Rev. 30, 310, 354, 408 (1927). 
® I. Langmuir, S. MacLane and K. B. Blodgett, Phys. Rev. 35, 478 (1930). 
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curves. The data for these curves were obtained by measuring J,, the electron 
current to the collector C, for low positive and negative potentials of C with 
respect to the point P. (The contact at the point L [See Fig. 1] had been so 
adjusted that it was at the same potential as the point P on the tungsten rib- 
bon for all values of the filament current.) The anode was maintained at a 
constant positive potential of six volts with respect to the point P during the 
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measurements. In Fig. 6 are shown Maxwellian distribution curves typical of 
the ones obtained in the experiment. The value of V at the point O was taken 
as the C.D.P. Theoretically the point O defines the state at which the point 
P and the collector are at the same potential.® Thus at this point the exter- 
nally applied potential is equal but opposite to the C.D.P. 


* N. B. Reynolds, Phys. Rev. 35, 163 (1930), 











88 HERBERT NELSON 


To determine the relation between V, and 7%, their values were deter- 
mined simultaneously for various values of @. 

The effect of field intensity on the C.D.P. and on the activity of the emit- 
ting surface was investigated by obtaining Maxwellian distribution curves 
for various anode potentials and by studying J,, the electron current to the 


collector, as a function of the anode potential V, at a constant potential of the 
collector. 


EXPERIMENTAL RESULTS 
I. The thermionic activity and the C.D.P. as functions of 0. 


The graph in Fig. 3 shows how V, varies with the time of evaporation. 
The temperature of the oxide coated filaments during the evaporation was 
estimated at 1200°K. 
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Fig. 3. Change in C.D.P, between collector and cathode with time of evaporation of barium. 





It is evident that there can exist a definite relation between the value of 
6 and the time of evaporation only when no migration'® of barium occurs on 
the tungsten surface. This condition appears to be fulfilled for @<1/2. For 
larger values of 0, migration seems to occur at the temperature at which the 
measurements are made. For this reason no attempts were made to obtain a 
quantitative relation between V, and 7, for values of @>1/2. On the other 
hand a qualitative relation between V, and 7» for the larger values of 6 could 
be obtained on the obvious assumption that 6 increases with the time of evap- 
oration. Thus it was found that V, increased and 75 decreased with @ until a 


10 Becker, “The Life History of Adsorbed Atoms and Ions” Proc. Am. Electrochem. Soc. 
55, 153 (1929). 
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maximum V, and a minimum 7’ were reached after which further evaporation 
caused V, to decrease and 7, to increase. This is in agreement with the re- ° 
sults obtained by Eglin" who found maximum activity of barium on tungsten 
for a value of @ he assumed to be equal to one, and lower activities for @>1. 

After the activity had passed the maximum, it reached a lower and prac- 
tically constant value which did not vary appreciably or regularly with 
further evaporation. This value of the activity was assumed to be that char- 
acteristic of barium metal. As the “thick” film of barium was removed from 
the tungsten ribbon by flashing it, the activity at first increased and then, 
after having reached a maximum, decreased as more and more barium was 
being removed until finally the state of clean tungsten was reached. 

The C.D.P. was always found to vary regularly with the activity and 
whenever the activity passed through a maximum the C.D.P. also passed 
through a maximum. 

Three separate determinations of the value of 7» for a “thick” film of ba- 
rium was made. The values obtained for the first, second, and third run were 
respectively 465°K, 465°K, and 460°K. The corresponding values obtained 
for V. were 2.0, 2.05, and 2.1 volts. For the maximum activity observed 

Te T 
/500\- (Tow) 





Gk) 


g 


w& 
r 


TEMPERATURE (DE 

















7) | 
q 4] / 2 3 
CONTACT POTENTIAL DIFFERENCE (VOLTS) 


Fig. 4. T is plotted as a function of the C.D.P. between collector and cat hode. 





the values for T, and V.. were found to be 410°K and 2.39 volts respectively. 
The values of 7, and V, for clean tungsten were 1450°K and —0.81 volts. 


II. Relation between V. and 7». 


The curve in Fig. 4 shows the relation between V, and 7 as obtained ex- 
perimentally. It was plotted from data obtained in three separate runs. The 
curve shows a straight line relation between V, and JT, between Tsy and 
Ton. Tow is the value of T» observed for clean tungsten and 79, is the lowest 
value observed for a “thick” layer of barium. At the point corresponding to 
To, there is a break in the curve, 7» varying less rapidly with V, for the values 
that are higher than for those which are lower than T9,. In Fig. 5 the region 

Eglin, Phys. Rev. 31, 1127 (1928). 
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of the break has been plotted on an enlarged scale. The points marked by cir- 
cles were obtained in the process of depositing barium while those marked by 
triangles were obtained in the process of removing barium. 

In the present experiment V, was measured to an accuracy of about 0.05 
volts. It is more difficult to estimate the accuracy of the temperature calcu- 
lations, especially in the lower range where corrections due to end losses must 
be made. In this region it is also possible that an appreciable error has been 
introduced in the temperature calculations due to the failure of taking into 
consideration the effect of the barium film on the radiation coefficient of 
tungsten. 
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Fig. 5. The region of the break of the curve shown in Fig. 4. 





In view of the above it follows that while the curve in Fig. 4 cannot be 
considered as giving the exact relation between V, and 74, especially at low 
values of 7» (high values of the activity), the observed linearity in the rela- 
tionship between V, and 7+, however, is sufficiently definite, and the break 
observed at Ton, sufficiently sharp to make the curve significant. 

It appears, on the other hand, that a more accurate investigation of the 
relation between V, and 7» simultaneously with determinations of the work 
function, ¢, and the A constant of the thermionic emission equation for va- 
rious values of 9, is desirable before attempts are made to interpret theoret- 
ically the results obtained. Such an investigation is now in progress in these 
laboratories. 


III. Velocity distribution determinations. 


The velocity distribution curves in Fig. 6 were obtained by plotting the 
logarithm of the current to the collector as a function of the potential differ- 
ence between the collector and the point P on the tungsten filament. Curve 
A is for a clean tungsten surface at temperature of 1423°K. Curve B was ob- 
tained for a barium on tungsten surface at a temperature of 700°K. 

On the assumption that the distribution of velocities of the emitted elec- 
trons is a Maxwellian one, Richardson shows that for parallel electrodes the 
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relation between retarding potentials and electron current may be expressed 
by the equation” 


i = ige Ve! AT 
where 7 is the value of 7 at V;=0. Vj is the retarding potential, e is the elec- 


tronic charge, k is the Boltzmann constant, and T is the temperature of the 
electron source. 
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Fig. 6. Typical Maxwellian distribution curves. 


The above equation may be written 
Ini = — Vye/kT + In wy. 


Thus if In z is plotted as a function of V1, the slope of the resulting straight 
line should be equal to e/kT. Since e and k are known constants it is evident 
that the temperature of the electron source may be determined from the slope 
of the velocity distribution curve. 

In Fig. 2 the temperatures corresponding to the points marked by crosses 
were calculated from Maxwellian distribution curves. These temperatures 
are on the average considerably higher than the true temperature of the fila- 
ment. Such a disparity between true temperature and temperature as calcu- 
lated from Maxwellian distribution curves has been a common observation in 
previous investigations of composite surfaces." 


IV. Effect of field intensity on activity. 


Curves such as are shown in Fig. 7 were obtained in studying the effect of 
field intensity on the thermionic activity for various values of 9. The data for 
these curves were obtained by measuring the electron current J, to the collec- 
tor for different values of the anode potential, the collector being maintained 


at a positive potential of 6 volts with respect to point P of the tungsten fila- 
ment. 


® Richardson, “Emission of Electricity from Hot Bodies,” p. 159. 
1 L. R. Koller, Phys. Rev. 22, 671 (1925). 
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The curve B was plotted from data obtained for a surface of approxi- 
mately the same activity as the one corresponding to curve A. The only dif- 
ference between the two cases is that the surface corresponding to curve A 
was obtained as a barium film was built up through evaporation from the 
oxide coated filaments, while the surface corresponding to curve B was ob- 
tained in the process of removing a “thick” layer of barium from the tungsten 
ribbon. 

It is evident that the variations in J, with the anode potential as exhibited 
in curve B cannot be due to change in activity. J, is, however, also affected 
by other factors which are functions of the anode potential. Secondary emis- 
sion from the collector and the aperture H are such factors. A third factor 
which appears necessary for the explanation of the great dissimilarity be- 
tween curves A and B is the presence in the electron current 7. of secondary 
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Fig. 7. The effect of the anode potential on the current to the collector at a 
fixed potential of the collector. 
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electrons which originate from points on the anode subjected to a compara- 
tively heavy bombardment of primary electrons from the partially coated 
tungsten filament and then are reflected through the aperture J7. Normally 
when the activity of the point P is greater than or nearly equal to that of the 
other parts of the tungsten filament, probably only a negligible fraction of J, 
consists of such electrons. When points exist on the filament, however, which 
are of considerably higher activity than the point P it is conceivable that a 
large fraction of J. consists of reflected secondary electrons. In such a case a 
decrease in J, with an increase in anode potential as exhibited in curve B may 
be accounted for as being due to the fact that for low anode potentials the 
secondary emission consists largely of reflected electrons of energy equal to 
that of the primary electrons while for higher anode potentials the secondary 
emission consists largely of low velocity electrons. Since the collector is main- 
tained at a positive potential of six volts it follows that secondary electrons 
which are emitted with velocities more than six volts lower than the velocity 
of the primary electrons will not reach the collector. Thus the fraction of the 
secondary electrons which have sufficient energy to reach the collector is 
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greater at lower anode potentials than at higher and 7. decreases as the po- 
tential is increased. 

Evidence of the presence of secondary electrons in J, as discussed above 
may be seen in the velocity distribution curves. Curve B in Fig. 8 corre- 
sponds to a case in which the activity of the point P on the tungsten filament 
was considerably less than that of a point directly below the oxide filament Q,. 
This condition had been obtained by glowing only the filament O; while or- 
dinarily both filaments were glowed to obtain a more uniform deposit of 
barium. The curve A corresponds to a case where the deposit was uniform. 
The gradual break of the curve B as compared to the sharper break of the 
curve A is attributed to the presence of secondary electrons. 

Non-uniformities in the barium film and resulting non-uniformities in the 
activity along the tungsten surface apparently were produced as a result of 
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Fig. 8. To show the effect of the presence of secondary electrons on the distribution curves. 


flashing the tungsten filament in removing the barium as well as a result 
of uneven deposit of barium in the process of evaporation. 

Due to the existence of the effects discussed above, the attempt to study 
the relation between field intensity and activity was abandoned. 

In the investigation of the relation between 75 and V, the assumption 
was made that the electron current J, had its source at the point P. In order 
to make that assumption valid the data were taken under conditions where 
the fraction of J, consisting of secondary electrons was negligible. A sharp 
break at the point O of the velocity distribution curve was taken as an indica- 
tion of a low fraction of secondary electrons. For most surface conditions 
studied, the fraction of secondary electrons was low. Conditions correspond- 
ing to the curve B in Fig. 7 and the curve B in Fig. 8 were rare. 

In conclusion, I wish to acknowledge helpful discussions during the course 
of this investigation with Dr. Saul Dushman and Dr. K. H. Kingdon of the 
General Electric Research Laboratory. I am also indebted to Dr. G. R. Shaw 
of this laboratory for helpful suggestions and advice. 


44 This explanation of the curve B is based upon suggestions offered by J. A. Becker. 
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THE EFFECT OF THE PRESENCE OF A GRID UPON 
CERTAIN CHARACTERISTICS OF THE AIRFLOW 
AT THE SURFACE OF AN AIRFOIL 


By Merit Scort* 
CoRNELL UNIVERSITY 
(Received May 18, 1931) 


ABSTRACT 


In a previous article it was reported that the heat dissipation per unit area from 
the surface of an airfoil, as measured by resistance strips mounted parallel to the span 
and flush with the surface, is characteristic of the airflow. An integral relation was set 
up from which it was shown that the coefficient of heat dissipation is proportional to 
the velocity gradient at the surface of the airfoil and hence to the viscous resistance. 
The effect of the presence of a grid upon the above mentioned characteristics has now 
been observed and comprises the subject of this paper. Experimental data are pre- 
sented to show the changes in the heat dissipation for different parts of the airflow due 
to grids of varying fineness and spacing, due to changes of angle of attack and due 
to changes in location of the grids with respect to the airfoil. In general, the heat dis- 
sipation is increased whenever a grid is present to induce turbulence in the air stream 
and the effect is maximum at about 0.6 x/c. It is concluded that near the leading edge 
there exists a region that serves rapidly to damp the turbulence. 


F Vpn to a development presented in an earlier paper! and as 
mentioned in the above abstract, the coefficient of heat dissipation is 
proportional to the velocity gradient at the surface of the airfoil and hence 
to the viscous resistance. This may be expressed as follows: 


Ou 00 
Mies afr aa} a) 
OVo Ovo 


Wherein: r=viscous resistance per unit area; u=coefficient of viscosity; 
0u/Oyo=velocity gradient .at surface; k = coefficient of thermal conductivity; 
00/dyo=temperature gradient at surface; and p.=coefficient of heat dis- 
sipation defined as the heat dissipated per unit area per degree difference in 
temperature between the air stream and the strip. The relation was shown to 
hold so long as the velocity gradient at the surface is not less than a certain 
limit, which might occur only in the immediate neighborhood of a stagnant 
point or a reversal of flow. 

The effect of the presence of a grid upon these characteristics has now been 
observed and comprises the subject of this paper. The measurements have 


been made upon a symmetrical Joukowski airfoil of 4’’ chord and 7” span, 


* Guggenheim Research Physicist. This work was performed in the Department of Physics 
at Cornell University under the direction of its committee on Aeronautic Research, supported 
by a grant from the Daniel Guggenheim Fund for the Promotion of Aeronautics. 

1M. Scott, Physics 1, 48 (1931). 
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with nine resistance strips mounted at 0.04, 0.155, 0.25, 0.36, 0.475, 0.575, 
0.69, 0.79, and 0.93 x/c. The airfoil was described and the heat dissipation 
characteristics presented in the former publication.! The mean working sec- 
tion velocity was maintained constant at 19.9 m/sec. 

The wind tunnel, in which the measurements were carried out, is dia- 
gramed in Fig. 1. The airfoil was mounted at W and the honey comb in ad- 
vance of the fan at 77. After considerable preliminary testing, it was found 
that satisfactory conditions for these experiments were obtained by omitting 
the section of the tunnel A BC and closing the throat as indicated in the figure. 
Measurements were made with this modification. There were no deflectors or 
straighteners between the entrance C and the airfoil at W, but the straight- 
ener at /H/ served satisfactorily to eliminate the angular swirling of the air 
stream. The span of the model was sufficiently less than that of the jet that 
the velocity variation across the span was not greater than 1.5 percent. 











Fig. 1. Diagram of the wind tunnel. 


The results of observations with several different grids are reported. 
These consisted of wires or screening extending across the entire jet and 
mounted at G(3-3/16’’ upstream from the leading edge) unless otherwise 
stated. The specifications of these grids are as follows: 


Grid No. 1 No. 14 copper wires (0.064’’) spacing 1/2’ 
Grid No. 2 No. 14 copper wires (0.064’’) spacing 1°’ 
Grid No. 3 No. 14 copper wires (0.064’’) single wire 
Grid No. 4 No. 26 copper wires (0.014’’) single wire 
Grid No. 5 1’’ mesh chicken wire 

Grid No. 6 No. 14 Copper wires (0.064’’) spacing 3/8’ 


In Fig. 2 are shown the results of measurements with grids No. 1, No. 2, 
No. 3 and No. 4. In all cases of single or parallel wire grids, the grids were 
mounted with the wires at right angles to the leading edge of the airfoil. In the 
case of single wires, these were mounted in the median plane of the airfoil, 
where the thermocouples for measurement of strip temperatures were at- 
tached to the strip.! At first thought one might be inclined to suppose that a 
grid would shield the airfoil and hence decrease the heat dissipation but this 
is clearly erroneous. On the contrary, the grid introduces turbulence, that is, 
a set of component velocities, oscillatory in character and superimposed upon 
those already present, characteristic of the tunnel and more particularly of 
that section of the tunnel at which measurements are made. This addi- 
tional turbulence increases the rate of transport of momentum to the surface 
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and the average velocity gradient at the surface. By Eq. (1) it is to be ex- 
pected therefore that the coefficient of heat dissipation should likewise be in- 
creased. 
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Fig. 2. The variation of the heat dissipation coefficient with percentage of chord abscissa 
both with and without four different grids. Mean working section velocity 19.9 m/sec. 

Fig. 3. The variation of the heat dissipation coefficient with percentage of chord abscissa, 
both with and without grid, the parameter being angle of attack. Mean working section velocity 
19.9 m/sec. @ with grid. O without grid. 

Fig. 4. The variation of the heat dissipation coefficient with percentage of chord abscissa, 
both with and without single wire grid, the parameter being position of grid as indicated by 
short, vertical lines. Mean working section velocity 19.9 m/sec. @ with grid. O without grid. 


Introducfion of turbulent velocity components into the hydrodynamical 
differential equations was made by Osborne Reynolds? in his attempt to 


2 O. Reynolds, Lon. Phil. Trans. (A) 186, 123 (1895). 
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determine the criterion for laminar flow. H. A. Lorentz,? employing the same 
method, shows that the hydrodynamical equations for the principal motion 
correspond to those for the actual or turbulent motion save that the stress 
components have additional terms due to the transport of turbulent momen- 
tum by the turbulent velocity components. The analogy to the transport of 
principal momentum by molecular velocity components is inviting. It is per- 
haps for this reason that in certain problems an apparent coefficient of vis- 
cosity has been introduced to include the turbulent effect and whose differ- 
ence from the actual coefficient of viscosity has been regarded as a measure 
of turbulence. 

Observation of the curves of Fig. 2 shows that not all parts of the airflow 
at the surface of the airfoil are equally sensitive to turbulence produced in the 
body of the fluid upstream. The largest effect, for these particular experi- 
mental conditions, is present at a value of x/c equal approximately to 0.6. 

The curves of Fig. 3 indicate the effect of the presence of a grid when the 
angle of attack is changed. The transitions were actually followed more closely 
than the representative curves of the figure show. Two features are to be 
noted. As the angle of attack is increased, the effect of the turbulence moves 
forward on the airfoil, as should be expected, in view of the fact that the stag- 
nant point likewise moves away from the leading edge along the face (under 
side) of the airfoil. Any point on the back(top) of the airfoil accordingly cor- 
responds, at higher angle of attack, to a point farther downstream so far as 
the flow is concerned. 

The second feature is that at a certain point, as indicated by the short, 
vertical lines on the figure, the turbulence ceases to have any further effect. 
Without doubt, this point corresponds to separation of the airflow from the 
airfoil surface and indicates that the turbulence travels downstream without 
penetrating the layer of separation. This point, as observed in this way, was 
found to move forward to about 0.1 x/c, then jump back to a point about 0.6 
x/c and gradually advance to approximately 0.4 x/c, as the angle of attack 
was increased from 0° to +25°. 

The next series of observations refers to introducing the turbulence actu- 
ally into the boundary layer next to the airfoil, using a single wire grid No. 3. 
The wire was located in the median plane of the airfoil, normal to its surface, 
and with the free end abutted to the airfoil surface. The results are given in 
Fig. 4, wherein the location of the grid is indicated by a vertical line neigh- 
boring the curve. In the bottom graph the grid was at the leading edge. There 
are several points to be noted in these data, which, for clarification, have 
been replotted in Fig. 5 in which the short, vertical lines indicates the leading 
edge. 

First, the increase in the heat dissipation coefficient due to the grid in 
general increases (see Fig. 5) as the grid approaches the strip, then either im- 
mediately damps to zero as for strips No. 1, No. 2, No. 3 and No. 4 or more 
gradually decreases to zero after the grid has passed the strip. Second, the 
maximum effect of the grid increases, as an examination of the figure shows, 


3H. A. Lorentz, Abh. u. theoret. Phys. p. 60, part I). 
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with percentage of chord to a value of the latter of 0.6, then decreases. This is 
likewise true of the fractional maximum increase of the coefficient. Third, it is 
to be noted that for strips No. 5, No. 6, No. 7, No. 8, and No. 9, the effect of 
the grid increases rapidly (see Fig. 5) after the grid has passed the leading 
edge and entered the boundary layer, but this increase continues only for 
about 0.4 x/c, when it has reached a maximum value. 
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Fig. 5. The variation of the increments of the heat dissipation coefficient with location of single 
wire grid with respect to strip for the various strips. Data taken from Fig. 4. 


The principal conclusion to be drawn is that the turbulence, induced by 
such grids, is rapidly damped in the boundary layer immediately downstream 
from the stagnant point (the leading edge in the case of angle of attack equal 
to zero). The damping decreases with distance back, until at about 0.4 x/c the 
effect seems to have vanished. The phenomenon is undoubtedly associated 
with and descriptive of the laminar flow in the boundary layer in the vicinity 
of the leading edge. 

The following experiments are now described to present a more intimate 
view of the influence of the grids. A single No. 14 copper wire is mounted 
parallel to the leading edge and at a point A 3-3/16’’ in advance of it, as 
shown in Fig. 6a. The increase in the average heat dissipation coefficient due 
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to it at x/c=0.575 (the thermocouple mounted at point C in the figure) is 
approximately 1.4 watts/cm?°C. If the wire is moved to point B, the effect 
of the grid becomes negligible. If, on the other hand, the same wire is 
mounted at right angles to the leading edge as at A, in Fig. 6b, the increase 
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Fig. 6. 


in the average coefficient is only 0.5 watts/cm?-°C but the wire must be dis- 
placed about twice as much, to B, before the effect of the grid becomes negli- 
gible. Grids of No. 14 copper wire of spacings 1’’, 1/2’’ and 3/8’’ are now 
interposed at the same point 3-3/16’’ in advance of the leading edge and ob- 
servations made with the wires parallel to the leading edge and perpendic- 
ular to it. Care must always be taken to keep the central wire in the plane of 
the airfoil and in the median plane respectively. The results are given in Table 
I wherein Ap, is the increase in the average heat dissipation coefficient due to 
the grid. 








TABLE I 
Abe watts/cm?-°C 
wires parallel to leading wires perpendicular to 
edge leading edge 
Single wire 1.40 0.50 
1” spacing 1.40 0.81 
1/2” spacing 1.40 1.46 


3/8” spacing 1.81 1.91 





The explanation of the increase in Ap, for the case of the grids with wires 
perpendicular lies in the fact that the wake of the wires does not cover the 
entire airfoil and the heat dissipation is increased only at those parts of the 
strips which lie in the wakes of the wires. The fact that Ap, does actually be- 
come somewhat larger for the perpendicular than for the parallel wires (and 
this has been rather carefully verified) is presumably due to the ability of the 
vortices with perpendicular axes to penetrate the boundary layer more freely 
than those of parallel axes. 

If wires of varying size are mounted perpendicularly, as at A in Fig. 6b, 
the result of measuring Ap,’s at x/c =0.575 are as shown in Fig. 7. The con- 
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sistency of the data for any single series without interrupting the air stream 
is to be noted, as shown by the dots or by the circles, whereas from one series 
to another, there seems to be a constant difference although the different 
parameters under control were maintained the same. This was characteristic 
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Fig. 7. The variation in the increments of the heat dissipation coefficient with radius of 
wire, located as at A, Fig. 7a and as at A, Fig. 7b. Dots @ and circles © represent single series 
of data without interruption of airstream, wires perpendicular. Circles with bars-O- represent 
single series with wires parallel. 


of the work throughout. If the wires are mounted parallelly, as at A in Fig. 6b, 
the corresponding results are as given also in Fig. 7 by the circles with bars. 

In conclusion, it may be remarked that the method herein described 
offers a sensitive means of observing turbulence in its relation to airfoil char- 
acteristics, and may be employed as a simple means of standardizing aero- 
dynamic tests relative to the parameter turbulence. 
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ABSTRACT 


The study of colloids has reached the point where many important contributions 
are coming from physicists. Even more attention from physicists will be required if this 
science is to progress as it should. A suitable dispersion of zinc oxide in kerosene has a 
high yield value and marked plastic properties. A few drops of a dispersing agent added 
to this mixture largely destroys its plastic properties and greatly increases its fluidity. 
It is well established that the dispersing agent accomplishes this result by decreasing 
the interfacial tension and thereby destroying the flocculation of the pigment. In cer- 
tain cases in paints, however, there is reason to believe an entirely different phenome- 
non is responsible for yield value and plasticity. Much work with the tools of the 
physicist is indicated in the solution of this problem. There is much to be done in the 
study of solid-liquid interfacial tension, adsorption layers, interfacial potentials, and 
the dielectric properties of colloids before these problems are to be solved. The be- 
havior of pigments in rubber is also discussed. 


{pe study of colloids has passed through the various stages which are 
common to most sciences. Once the general principles underlying colloid 
phenomena were discovered, the more salient features of a great number of 
different colloid systems were studied and described. Empirical rules for pre- 
paring, studying, and describing these systems were evolved, and the quantity 
of more or less disconnected data which came to exist was quite overwhelm- 
ing. Then, gradually, came the tendency to return to the more simple systems 
and study them exhaustively. With this reaction came new methods of ex- 
perimental attack and new viewpoints on the most fundamental and estab- 
lished principles of the science. As these expermental methods and these new 
viewpoints became more and more fundamental, they tended to change, as 
always, from those traditional to the chemist to those traditional to the 
physicist. Many colloid chemists became physicists and some physicists be- 
came interested in colloids. The latter was probably the more efficient transi- 
tion; it is easier to apply familiar methods to unfamiliar problems than un- 
familiar methods to familiar problems. The contributions of certain physicists 
to the field of colloids have been notable. 

Pigments are used in paint and rubber largely to influence the chemical 
and mechanical properties of the systems. The problems arising in the study 
of these disperse systems are probably more numerous and complex than 
would be expected by the uninitiate. While much progress has been made in 
the manufacture and application of these pigments, many of the most funda- 
mental properties of the dispersions are not at all well understood. The 
writer wishes to present a few of these problems to the group which this pub- 
lication will reach, with the hope that perhaps a few physicists will become 
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sufficiently interested to apply their talents to these most interesting and im- 
portant questions. 

When a mixture of a nonpolar oil, such as kerosene, and zinc oxide is 
prepared in suitable proportions, a highly plastic material is formed. This 
paste has a relatively high yield value and is capable of very little elastic 
deformation; under higher stresses it flows continuously. Such properties 
would hardly be conceivable if every particle were completely surrounded 
with liquid; as a matter of fact, it is not. The free energy of the interface is so 
high that a decrease in free energy results when the solid particles approach 
closely enough to one another to eliminate some of the interface. In a dilute 
suspension this process would result in a clumping of the pigment; that is, 
the pigment would be flocculated. In this concentrated paste the pigment 
forms an essentially continuous structure in which distortion involves what 
amounts to solid-to-solid friction. Thus is the yield value and the absence of 
marked elasticity explained. 

If a few drops of a surface active material, such as stearic acid, be added 
to the mixture, it will concentrate at the interface, lower the interfacial ten- 
sion, and reduce or eliminate the tendency of the pigment to flocculate. Thus 
each particle becomes entirely surrounded with a true liquid layer and the 
yield value disappears or at least is reduced to a different order of magnitude. 
These phenomena are well-known and the theoretical analysis of them is well 
substantiated and quite generally accepted. There are, however, other pig- 
ment-liquid mixtures with somewhat similar properties, which have received 
no such satisfactory theoretical treatment. 

A suitable ferric hydroxide-water gel, which has a considerable yield 
value, may be transformed to a sol by shaking and changed back to a gel by 
by being allowed to stand. This phenomenon is called thixotropy. Such a 
system differs in properties from the kerosene-zinc oxide mixture in two im- 
portant respects. First, if stress in excess of the yield value has been applied 
to the ferric hydroxide sol, a sensible length of time is required for the yield 
value to be regained while in the case of the zinc oxide-kerosene mixture the 
regain is essentially instantaneous; and, second, the ferric hydroxide sol is 
capable of a marked elastic deformation while the other system is not. 
While these differences may be a matter of degree rather than kind, it is 
evident that the theory which was so satisfactory for one system can hardly 
be applied to the other. A system whose yield value is the result of essentially 
solid to solid friction could not give the large elastic deformation under a 
stress below the yield value that is characteristic of a gel; and it seems 
hardly probable that much time could be consumed in assuming a stable 
flocculated structure in a medium so fluid as water when it is considered that 
the change in the position of each individual particle in going from a com- 
pletely dispersed to a flocculated state would be only a micron or so. And, to 
add to the confusion, even though a satisfactory theoretical treatment of 
the gel were available, there are systems with border-line properties which 
could not, at present, be ascribed definitely to either class. 

One hypothesis has been widely advanced to cover elasticity and yield 
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value in thixotropic gels. In a word, this hypothesis supposes the existence 
of an adsorbed layer of dispersing medium around each solid particle so thick 
that relative motion of the particles is impossible without at least partial 
destruction of the adsorbed layer. The elasticity inherent in this envelope 
accounts for the elasticity of the gel and the time required for the formation 
of the layer accounts for the duration of the sol state. This picture fits the 
facts well enough to have gained considerable cognizance. The difficulty lies 
in the fact that the existence of this envelope has never been conclusively 
demonstrated; has never, in fact, so far as the writer is aware, been shown to 
be compatible with what is known about the molecular forces and the kinetics 
involved. 

To account for the data that have been published by some workers! the 
adsorbed layer of water around the micelles in certain gels would have to be 
some microns thick. If these water molecules are fixed firmly enough to form 
a structure of the observed mechanical strength, the phenomenon should be 
accompaied by certain other effects. In the first place, the constraint of so 
many molecules would be expected to be accompanied by marked heat effects, 
and certainly the heat capacity of the constrained water should be lower 
than that of the free liquid. Some data on these effects would be very helpful 
at the present time. 

It has also been pointed out that polar molecules so constrained should 
respond less readily to an electrostatic field than those of the free liquid, thus 
reducing the dielectric constant of the material. Since the molecules in the 
adsorbed layer are not in a truly solid condition but more likely simply have a 
much increased relaxation time, this effect should be especially pronounced 
at high frequencies. A number of attempts have been made to determine the 
bound water in certain systems but most of them have been concerned with 
systems like gelatine which are a bit beside the present discussion. An excellent 
paper by Kistler? has appeared very recently, however, on the dielectric 
constants of Fe.O3;, Al,O3;, and bentonite sols. These sols are all non-thixo- 
tropic; that is, they will not gel if they are sufficiently free of electrolyte. Now 
Kistler made the very surprising observation that these sols had dielectric 
constants indistinguishable from that of water in the non-thixotropic state 
but that the di-electric constant was lowered when they were rendered thixo- 
tropic by the addition of electrolyte, whether the measurement was made on the 
sol or the gel. The decrease in the dielectric constant was a function of the 
amount of electrolyte added but was constant for a given material whether 
it was measured as a gel or a sol. This is a most striking observation and it 
is very desirable that this work be checked by some other worker and ex- 
tended to other systems. 

The writer, in collaboration with Mr. S. E. Kamerling, made a sporadic at- 
tempt to obtain some dielectric constant data on thixotropic and non-thixo- 
tropic paints nearly a year ago. The experimental difficulties were hardly 
overcome in the time available and no conclusive data were obtained. The 

1 E. A. Hauser, Kolloid-Z. 48, 57 (1929). 

2S. S. Kistler, J. Phys. Chem. 35, 815 (1931). 
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study of non-aqueous systems offers certain advantages as well as certain 
difficulties that are not found in water dispersions. Any well substantiated 
data in this field would certainly be a valuable contribution. 

Early in this paper a statement was made in regard to the relationship 
between flocculation and interfacial tension. In many cases it is difficult to 
determine whether the tendency to flocculate has been increased or decreased 
by an alteration in a disperse system which has resulted in a change in its 
physical properties. For instance, certain small changes in the composition 
of a pigment-vehicle mixture may greatly alter its flow characteristics under 
very small stresses. Whether this is a result of a change in interfacial tension 
with a consequent change in the flocculating tendency or to a change in the 
thickness of the adsorbed layer, is frequently quite impossible to establish. A 
method whereby solid-liquid interfacial tensions could be determined would 
be of enormous value here as in other phases of surface chemistry. Failing 
this, the determination of differences in interfacial tension from one system 
to another would be a big step. It is recognized, of course, that the difficulty 
in determining the tension of the solid-liquid interface is much the same as 
that of determining the surface energy of the solid, a problem which has 
received the attention of a number of physicists. There are, however, a 
number of possibilities which may in time lead to some valuable data on this 
problem. 

Bartell and Osterhof* have designed an apparatus whereby it is claimed 
that adhesion tension measurements can be determined. This makes it possi- 
ble, in certain cases, to determine the difference in the interfacial tension 
between two liquids and the same solid. There has been some question as to 
the validity of the data obtained by this method, however, and the difficulty 
has not been entirely cleared up. 

There are two other possible attacks which might shed some light on this 
problem. Koch‘ has derived the formula, 


20M 
ronk 


Ae = 





where Ae is the difference in potential between a spherical particle of radius 
r and a plane surface of the same material and @ is the interfacial tension 
between the material in question and the liquid in which it is immersed. Since 
the other variables in the equation are capable of experimental determination, 
the surface tension could be calculated. Although this expression is closely 
related to the change of solubility with interfacial tension, the experimental 
approach would be quite different. The practicability of this method would 
be well worth investigating. 

In many systems flocculation and dispersion are reversible phenomena; 
that is, a flocculated system will disperse spontaneously if the interfacial 
tension is sufficiently reduced, and vice versa. Where the interfacial tension 
is so low that no particles adhere upon collision, the dispersion is complete. 


3 Bartell and Osterhof, Ind. & Eng. Chem. 19, 1277 (1927). 
4 F. K. V. Koch, Phil. Mag. Sup. 2, 585 (1931). 
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When, however, the interfacial tension is somewhat larger, the tendency to 
flocculate due to the adhesion of the particles must reach an equilibrium 
with the tendency of the particles to disperse due to their thermal agitation. 
It should be possible to formulate the energy relationships involved at equilib- 
rium and, by experimentally determining the degree of dispersion, to calcu- 
late the interfacial tension. To the writer’s knowledge, no such approach to 
the problem has been made. 

There are many problems in connection with the dispersion of pigments 
in rubber that might well receive the attention of physicists, but the field 
can only be touched upon here. The addition of pigments, such as zinc oxide 
and carbon black, to rubber causes a marked reinforcement in the cured 
stock; that is, a pigmented stock has a stiffer stress-strain curve and requires 
a greater total energy to break than an unpigmented stock. A great deal of 
work has been done in this field, but the facts still remain without a satis- 
factory explanation. The separate effects of particle size, particle shape, and 
the specific properties of the material composing the pigment have never been 
adequately treated. This situation is probably the result of poor fundamental 
approach as well as the inherent experimental difficulties attending quantita- 
tive studies of rubber compounds. 

One example of faulty approach to this problem may be worth mention- 
ing, largely because the theory which it involves seems to be receiving more 
general acceptance than the facts justify. The theory is based upon the follow- 
ing statement by Hock,® ‘“‘Man hat gefunden, dass diese verstiirkende Wir- 
kung der aktiven Fiillstoffe, zu denen z.B. Russ and Zinkoxyd gehéren, um 
so auffilliger hervortritt, je feiner die Fiillstoffe sind. Die Benetzung von 
Kautschuk und Fiillstoff is also offenbar die Ursache der Verstiirkung. Sucht 
man auf diesen Grundvorstellungen eine Theorie der Kautschukfiillstoffe zu 
entwickeln, so hat man offenbar die freie Grenzflichenenergie zwischen Kaut- 
schuk und Fiillstoff fiir die Vergrisserung der Zerreissarbeit verantwortlich 
zu machen.” From this standpoint is the theory developed. Hartner,® having 
devised a method for estimating the total free energy of the pigment-rubber 
interface, gets an order-of-magnitude check with the increased energy to 
break. The point which both of these workers seem to have overlooked is that 
this free surface energy can be of primary importance in reinforcement, only 
if the interface changes in magnitude or character during elongation and 
break. It has been established experimentally many times that the former 
change does not take place, and there is no reason to believe that the latter 
does. As a matter of fact, Hartner seems to assume that the entire rubber- 
pigment interface is destroyed when the test piece is ruptured, an obviously 
erroneous assumption. This case is cited as an example of the many practical 
problems concerning the properties of pigment dispersions which would 
doubtless progress more rapidly if they could receive the benefit of the ex- 
perimental technique and method of attack of the physicist. 


5 L. Hock, Z. Elektrochem. 34, 662 (1928). 
6 F, Hartner, Kolloidchem. Beihefte 30, 83 (1929). 
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ABSTRACT 


This report refers in particular to suspensions of biological cells. The interior of a 
biological cell is a well-conducting salt solution, but at its surface polarization or the 
presence of a nonconducting membrane obstructs the passage of an electric current 
with a resulting low conductivity and high dielectric constant of the suspension. The 
theory for the passage of an alternating current through a suspension of this type is 
dealt with and experimental work with various kinds of suspensions such as cream, 
wet sand, blood, and animal tissue, is described to show the applicability of the theory. 
The determination of the volume concentration of suspensions, of the form factor of 

¢ the suspended particles, of the thickness of surface membranes, and of the interior con- 
ductivity of living cells, are instances of practical applications. 


T HIS report includes a group of investigations which have been made dur- 
ing the last few years.! The work was done with a special interest in 
biological usage, but applies as well to a wider field. 

The first measurements of the electric resistance of living cells were made 
in the middle of the foregoing century. Such well-known physicists as Peltier, 
Weber, v. Humboldt, and Lenz worked in this field. A large amount of 
work was done by the noted physiologist du Bois Reymond, who summarized 
the work done up to 1850 in his “Untersuchungen iiber tierische Elektrizitat.” 

The fundamental fact was brought out that a living cell becomes polarized 
when an electric current passes through it, the resistance of living matter 
being, in consequence, high when a direct current or a current of low fre- 
quency is employed, but decreasing when the frequency is increased. Very 
little of lasting value was accomplished, however, before the advent of the 
audion oscillator made the production of alternating currents of pure form 
and of different frequencies a simple matter. 

For theoretical study, we may take as a simplified model a suspension of 
conducting particles, at the surface of each of which polarization of the elec- 
tric current takes place. Consider first the conductivity of a suspension of 
homogeneous, non-polarizable spheres. The theoretical solution of this prob- 
lem was given by Maxwell in his “Treatise on Electricity and Magnetism.” 
He derived the well-known formula: 

(K/K,)-1  (Kx/K,) —1 es 
(K/K,) +2” (K2/K,) +2 
1 Phys. Rev. 24, 575 (1924); 25, 361, (1925); 26, 678, (1925); 26, 682, (1925). 
Jour. Gen. Physiol. 6, 741 (1924); 9, 137 (1925); 9, 153 (1925). 
Jour. Cancer Res. 10, 340 (1926). 
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where K =Specific conductivity of suspension; K,=Specific conductivity of 
suspending liquid; K»=Specific conductivity of suspended particle; and p 
= Volume concentration of suspended material, which in the field of conduc- 
tivity is usually referred to as Maxwell’s formula, but which is known under 
other names in other departments of physics, such as the Lorenz-Lorentz for- 
mula, the Clausius-Mossotti formula, or the Poisson formula. 

I have tested this formula over a wide range of volume concentrations 
by working with cream. A heavy cream was obtained by separating raw, 
fresh milk in a De Laval cream separator. The conductivity of this cream, of 


TABLE I. Test of Maxwell's formula for the conductivity of a suspension of spheres. 


Volume concentration percent 





Calc. from Calc. from 

resistances dilution factors 
Primary cream 62.73 +.12 62.93 + .19 
Dilution No. 1 40.40 + .08 40 .44+.12 
Dilution No. 2 22.49 + .06 Standard 





the skimmed milk, and of various mixtures of the two were measured. Table 
I shows the result of such a series. One column gives the volume percent of 
the fat particles calculated from Maxwell’s formula, which in this case reads: 
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Fig. 1. Resistance of suspensions of red corpuscles of different concentrations. 


where K=Conductivity of cream; K,= Conductivity of skimmed milk; and 
p= Fat-content of the cream. The next column shows the relative volume con- 
centration calculated from the dilution factor. 

Much work has been done on the conductivity of blood. The red corpuscle 
is a non-conductor and therefore the conductivity of blood is lower than the 
conductivity of the serum. 

How well does Maxwell’s formula apply to a suspension of red corpuscles? 
In Fig. 1 are plotted experimental data obtained with red corpuscles of a dog. 
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The ordinates indicate volume concentration, and the abscissas the ratio of 
conductivities of serum to corpuscle-suspension. The upper curve is calcu- 
lated from Maxwell’s formula. The reason for the deviation is that the form 
of the corpuscle of mammals is far removed from the spherical shape, it 
being that of a flat, biconcave disk. 
It is a fairly simple matter to extend Maxwell’s formula to the case of a 
suspension of spheroids. The following formula is obtained: 
(K/K,) — 1 (K2/K,) —‘1 ( 
Se 3) 
(K/K,) + Xa (Ko/K,i) + X 
where X depends on K2/K, and on the axis ratio a/b of the spheroid, but is 
independent of the volume concentration p. X is shown graphically in Figs. 
2 and 3 for the oblate spheroid and the prolate spheroid, respectively. When 
a/b=1, we have X =2, which is the Maxwell formula. 
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Fig. 2. Graphical representation of x for the case of the oblate spheroid. 


The curve II in Fig. 1 has been calculated from formula (3), using a/b 
=1/4.25, which gives X = +1.05. The agreement is very satisfactory; in 
fact, it is rather surprising that the agreement is so good at the highest volume 
concentrations, where the corpuscles are packed far beyond the limit set by 
the total packing factor, indicating that deformation of the cells must have 
taken place. There has been much discussion about what the form of a red 
corpuscle is when it is in the blood. It is, of course, not certain that the form 
we observe microscopically is the true form, and it has been insisted, from 
various sides, that it is not. In the microscope the corpuscle is seen to be of 
the form of a biconcave disk. The ratio of diameter to maximal thickness for 
the dog is between 3 and 5. Our value a/b =1/4.25 indicates this to be the 
true form. 

As a precision method for determining the red corpuscle volume, the 
conductivity method is not wholly satisfactory for mammals, because the 
form of the corpuscle is not an ellipsoid. But for the lower orders—birds, fish, 
reptiles—the shape is ellipsoidal and the conductivity method should give 
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results of highest exactness. The method is often used in physiological inves- 
tigations and might be useful also in practical clinical work if a simple appara- 
tus for this particular use were placed on the market. The conductivity 
method may also be useful for the practical determination of the fat content 
of milk and cream, although no attempt of so using it has been made as yet. 
Other possible fields of practical application are the determination of rubber 
in latex and the analysis of soils. 
TABLE II. Conductivity of suspensions of sand. 


Db. / : p (obs.) x 
ki/k (obs.) (percent) (calculated) 
; 1.190 10 1.410 
1.426 20 1.420 
1.734 30 1.405 
2.156 40 1.367 
2.715 50 1.400 
3.613 60 1 do 
Average value 1.392 


Table II shows some measurements of the conductivity of suspensions of 
sand, which were made in Wilhelm Ostwald’s laboratory many years ago. The 
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Fig. 3. Graphical representation of x for the case of the prolate spheroid. 


suspensions were made by adding the sand to hot salt solutions containing 
3 percent gelatine and cooling this mixture to gelatination under continuous 
rotation. The concentrations vary from 10 percent to 60 percent. I have ap- 
plied formula (3) to these measurements, for each concentration calculating 
X, using the experimental values of Ki/K and p. The values of X stay per- 
fectly constant. As X is smaller than 2, the sand particles must have been 
flat. The average ratio of a/b is about 1/3. 

For a current of low frequency the living cell is practically a nonconductor, 
due to polarization at its surface. In order to obtain the true conductivity of 
the cell, a current of high frequency must be used. This is illustrated by data 
shown in Table III. Here are given measurements of the resistance of calf 
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Tase IIT. Capacity and resistance of blood of calf for frequencies from 87,000 to 43 million cycles. 


146 1 (146 -10-" - w- 370)? 
Cycles per 370)? 4 


f Capacity Resistance 1+(146-107 370 1+(146-10-"-w-370)? 
rae (mmf) (ohms) 1 —1 
(w/27). 
191 
87,000 146 191 146 191 
833,000 130 181 135 183 
1.17-10 118 174 126 178 
1.52-— 106 168 115 172 
2.04-— 90 159 98 163 
3.04-— 68 148 71 151 
3.82-— 60 144 55 144 
4.52-— 39 138 43 140 
(124) 126 
re — extrapolated 


blood with currents of different frequencies. The resistance at infinite fre- 
quency is extrapolated to be 124 ohms, while at low frequency the resistance 
is 191 ohms. Using formula (3) with a/b =1/4.25, the value of 124 ohms leads 
to a value for the specific resistance of the corpuscle as follows: 
spec. res. corpuscle 
——— ee Ow SY 
spec. res. serum 


It may be noted that the high resistance of the red corpuscle, which is re- 
corded when the frequency is low, seems to be due not to a polarization at its 
surface but to the presence of a thin nonconducting membrane around the 
corpuscle. The interior of the corpuscle is a well-conducting salt solution. Ata 
sufficiently high frequency the impedance of the static capacity of the mem- 
brane becomes negligible; the resistance of the red corpuscle calculated above 
representing the resistance of the interior of the corpuscle. By this method a 
continuous record may be made of the exchange of electrolytes between the 
serum and the corpuscle interior,—a problem of considerable physiological 
importance. 

The polarization which takes place when an electric current is sent 
through living cells makes itself felt when we try to measure the resistance 
of the cells in a bridge. It is impossible to obtain an equilibrium without 
adding a certain amount of capacity in parallel to the resistance in the oppo- 
site arm of the bridge. For small current densities, a suspension of living cells 
is electrically equivalent to a resistance in parallel to a capacity |Fig. 4(c) |. 
The value of this parallel capacity we term the capacity of the suspension. 
The capacity of a cm cube of the suspension is the specific capacity. The speci- 
fic capacity of animal or vegetable tissues, such as liver, muscle, roots or fruits, 
is usually quite large—several thousand micromicrofarads at 1000 cycles— 
while the specific resistance is a few thousand ohms. The measurement of the 
capacity for such material is very simple. Suspensions of single cells, such as 


blood or bacteria, usually have much lower capacities. Since their resistance 
is low, the exact measurement of the capacity is a more difficult matter. For 
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instance, the specific capacity of the blood of man, which is due to the red 
corpuscles, is a few hundred micromicrofarads, while the specific resistance is 
about 150 ohms. Table 1V shows measurements at different frequencies of the 


TABLE IV. Resistance and capacity of suspension of red corpuscles of chicken. 











| Directly | Capacity | Capacity due Capacity 
Kilocycles | measured measured for to polarization| of corpuscles 
per | Resistance | capacity | salt solution | at electrodes | 
second | | (Micromicro- | (Micromicro- | (Micromicro- | (Micromicro- 
| farad) | farad) farad) farad) 
1/2 | 4080 | 366 | +46 56 321 
1 4080 340 +12 20 328 
2 | 4080 317.5 | —4 | 7 320.5 
4 4080 307.5 | 14.5 | 2.5 322 
8 4073 298.5 —17 1 385.5 
16 4022 B93 5 — 20 0 533.S 
32 3924 “293.5 —16 0 309.5 
64 3664 274 —16.5 | 0 290.5 
128 3072 231 —16 0 247 
256 2324 173.5 —16 0 189.5 
Temperature: 18.30°C Specific capacity: 825 micromicrofarads. 
Cell constant: 2.58 Specific resistance of serum: 92.8 ohms. 
Specific resistance: 1582 ohms Specific capacity for 100 percent suspension: 878 micromi- 


crofarads. 


capacity of a concentrated suspension of red corpuscles of chicken made with 
a Wheatstone bridge. The first column gives the resistance, the next the capac- 
ity directly recorded in the bridge. The next column gives the capacity meas- 
ured when the electrolytic cell was filled with a salt solution of the same 
resistance as the corpuscle suspension. The difference is the capacity due to 
the red corpuscles. The capacity of the salt solution increases at low frequen- 
cies due to electrode polarization. I have chosen this example to show the 
influence of electrode polarization. However, for accurate results we should 
use electrolytic cells so large that the electrode polarization is negligible. 

Consider a suspension of conducting spheroids, each being covered with a 
thin nonconducting membrane. The specific capacity of the suspension at low 
frequency is expressed by the following formula: 


’ K ‘ 


where K = Conductivity of the suspension; K, = Conductivity of the suspend- 
ing medium; g= Major axis of the spheroid; Co=Static capacity of one sq. 
cm of the membrane; and a= Depends on the axis ratio of the spheroid, and 
its values are shown in Table V. 


TABLE V. Values of a as a function of b/a and of a/b. 


b/a 1 2 3 4 


oo 
a 1.50 1.30 $27 1.28 1.65 
a/b 1 2 3 4 oo 
a 1.50 1 .( 


)3 0.94 0.94 (0.118 Xa/b) 


Axis of spheroid a, b, }, 
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We may write formula (4) as follows: 


K 
C= Cumti-( 1 -_ =) 


With Cio0% = agCo. 


(5) 


where Cioox, is a constant, which may be considered as the capacity of a 100 
percent suspension. This formula shows that the capacity varies as 1—K/K,, 
in its dependence on the volume concentration. In this form, the formula 
is independent of the form of the spheroid and may be applied to a suspension 
of particles of any form. 

This formula has been tested by measuring the capacity of suspensions of 
red corpuscles of different concentrations. Table VI shows the results of a 


TABLE VI. Capacity of suspensions of red corpuscles of dog in own serum, 





























Volume | Capacity 
| concentra- | (Cioo) for 
No. of tion cal- | Resistance; Capacity |100 percent Origi 
experiment! culated | (r) (ohms) | (C) mmf | concentra- — 
from tion (cal- 
resistance | | culated) 
percent 
I 83.9 931. 374 411 Concentrated by centrifugation 
from original blood. 
II 21.0 126.9 | 129 385 From 83.9 percent suspension by 
dilution. 
Ill 72.0 498. 343 411 From 83.9 percent and 21. percent 
' suspensions. 
IV 47.5 230.2 237 374 From 72.0 percent suspension by 
dilution. 
V 60.2 329. 286 385 From 83.9 percent and 47.5 per- 
cent suspensions, 














series of such measurements. By formula (5) Ciooe, is calculated from the 
observed values of the capacity C. The values of Cioo% are constant from 21 
percent to 84 percent concentration. 

Using formula (4) we find the capacity per cm?* of the corpuscle surface 
Co=0.80uf, the value obtained being the same for corpuscles of man, rabbit, 
calf, dog, chicken and turtle. 

An interesting fact about the red corpuscle is that its capacity is indepen- 
dent of the frequency of the electric current. The capacity also stays constant 
when the corpuscle is suspended in different media which do not change its 
form, such as a sodium chloride solution or a sugar solution. These facts 
favor the assumption that the capacity represents an actual static capacity 
of a thin non-conducting surface membrane. If we assume this to be true and 
choose a suitable value for the dielectric constant of this membrane we may 
estimate its thickness. With Cp) =0.80uf and the dielectric constant 3, we find 
the thickness to be 40.10-§ cm. 

Formula (4) holds only when the frequency is so low that the impedance 
of the capacity of the membrane is large compared with the resistance of the 
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cell interior. When the frequency is increased to such a point that this is not 
the case, the resistance and capacity of the suspension decrease. With good 
approximation this decrease may be calculated as follows: Fig. 4 (a) shows 
the suspension schematically. A part of the electric current passes between 
the cells and encounters the resistance Ro. A part passes through the cells and 


y C5 ct 


R 
Ro R; Ry 


(a) (b) (c) 


Fig. 4. Electrical diagram of a suspension. 


encounters the capacity of the membrane C, and the resistance R;, which is 
the resistance of the cell interior. Thus, electrically, the suspension may be 
represented by the diagram (b). The capacity and resistance of the suspension 
are defined by diagram (c). The following equations are easily derived: 


Co 
C= (6) 
1 + (CywwR;)? 


1 1 1 (CwR,)? 


R Ro Ri 1 + (CwR,)? 








(7) 


Co, Ro, and R; are all directly observable: Ro is the resistance measured at low 
frequency; Co is the capacity measured at low frequency; (1/Ro+1/R;)"! 
is the resistance at very high frequency. 

Table III contains at different frequencies measurements of the resistance 
and capacity of a suspension of red corpuscle, and the values calculated 
by formulae (6) and (7). The agreement is quite satisfactory. 

Usually the capacity of living cells is due to a polarization at the cell 
surface. The polarization acts as a capacity in series with a resistance and 
the resistance and capacity of the suspension may be expressed by formulae 
(6) and (7) by using: 


R; = Rp +R? 
Co = Co 


where R, is the polarization resistance, R;1 the resistance of the cell interior, 
and C, the polarization capacity. Thus the following formulae are obtained: 








1 1 C,wR; + m ; 
— a 1 po (8) 
R Ro 1+ (C,wR; + m)? 
C 
C= : (9) 
1 + (C,wR; + m)? 
m = CwR,. 
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For low frequencies when R; is negligible compared with the impedance of 
the polarizations capacity, the equations may be written: 





1 1 
— = — + Cwm (10) 
R- Ro 
& 
C =——— (11) 
1 + m? 


Experience shows that m is independent of the frequency, or at least approxi- 
mately so, a law which is known to hold also for the polarization at metal 
electrodes in aqueous solutions. 
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Fig. 5. Resistance and capacity of muscle of rabbit for currents of different frequencies. 


The points in Fig. 5 represent measurements of the resistance and capac- 
ity of a piece of rabbit muscle. The decrease of the resistance at low frequen- 
cies is due solely to the polarizations resistance, and may be accounted for by 
formula (10). Hereby m may be calculated. We find m =0.46. By formula (11) 
C, may now be calculated, and thereafter by using formula (4), the value of 
the polarization capacity per cm? of cell surface may be calculated. 

It has been found for the polarization at the so-called polarizable metal 
electrodes that the polarization capacity varies nearly exponentially with 
the frequency C,~w~*, and this same law is found to hold also for polariza- 
tion at cell surfaces. The exponent falls between 0 and 1/2 and is usually about 
0.30. In this particular case it is 0.28. 

The broken curves in Fig. 5 have been calculated by formulae (10) and 
(11), without considering the resistance of the cell interior, using m =0.46 and 
taking C, to vary as w~®*8, The resistance of the cell interior makes its 
influence felt at about 5000 cycles. The fully drawn curves have been calcu- 
lated with R;=250 ohms. The agreement is very satisfactory. 

The measurement of the polarization of living cells is one of the important 
physical methods for obtaining information about the nature of cell mem- 
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branes. It should also be noted that the existence of polarization forms a 
criterion for the presence of living cells, which in some cases may be of use 
when optical methods fail. It is a significant fact that the polarization seems 
to be of the same type as is observed at so-called polarizable metal electrodes in 
aqueous solutions, e.g. platinum in sulphuric acid. Since it seems certain that 
definite potentials are maintained across cell membranes, one might have ex- 
pected to observe the type of polarization which is produced at the so-called 
non-polarizable metal electrodes as silver in silver nitrate and which is char- 
acterized thereby that the polarization capacity varies inversely as the square 
root of the frequency. But apparently the mechanism by which the bioelec- 
trical potential is maintained is too slow to be effective within the period of 
one cycle of the alternating current. 
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ABSTRACT 


Three different physical methods for studying the capillary structure of porous 
materials are developed and applied to the study of the fine effective continuous capil- 
lary structure of wood. The three characteristic capillary properties measured are, (1) 
the ratio of effective capillary cross section to the effective length, (2) the average 
effective capillary radius, and (3) the maximum effective capillary radius. Some of the 
important characteristics of the porous structure of wood are briefly presented in order 
to show that the major part of the void structure, the fiber cavities, have very little to 
do with the permeability of wood. The pores in the pit membranes which connect the 
adjacent fiber cavities are shown to be the structure effective in controlling permeabi- 
lity of wood. The pores in the pit membranes in general, are below microscopic visi- 
bility, thus necessitating other means of study. The ratio of the effective capillary cross 
section to the effective capillary length is determined from electrical resistance meas- 
urements of strong salt solutions filling the void structure and from the specific resis- 
tance of the salt solutions in bulk. Strong salt solutions are used in order to make the 
surface conductivity negligible as well as the conductivity of the swollen cell wall. The 
average effective capillary radius is determined by combining the results of the pre- 
ceeding measurements with results obtained from hydrostatic flow studies of a con- 
tinuous water phase extending through the section. Poiseuille’s law is used for the 
calculations. A differential pressure drop apparatus was devised for these measure- 
ments in which the flow of liquid through a standard capillary tube and the test 
specimen connected in series is determined by measuring the pressure drop through 
each. The maximum effective capillary radius is determined by measuring the gas 
pressure required to overcome the effect of the surface tension of water in the capillary 
systems. The data show the large differences between the fine capillary structure of 
heartwood and of sapwood. 


4 pr study of the capillary structure of porous materials is not only of 
theoretical interest but also of considerable practical importance. For ex- 
ample, a more nearly complete knowledge of the fine capillary structure of 
wood should lead to better seasoning practice, improved impregnation with 
preservatives, and improved penetration of chips in chemical pulping pro- 
cesses. 

, Variations in the coarser capillary structure of wood, notably the fiber 
cavities, the pores, and the resin ducts, which make up the major part of the 
porosity of wood, have been shown to effect the permeability but slightly. 
On the other hand, the permeability of two specimens that have a coarse 
capillary structure practically identical may differ tremendously. (Table II). 


1 Maintained at Madison, Wis., in cooperation with the University of Wisconsin. 
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This is due to the fact that the fiber cavities are closed at both ends, the only 
communication from fiber cavity to fiber cavity being through the membranes 
of the bordered pits or through the cell wall. It is the size and also the number 
of the finer openings that control the permeability of wood. 

Fig. 1 shows a photomicrographic section of a typical softwood cut across 
the fibers. A bordered pit connecting two adjacent fiber cavities appears in 
cross section below the center of the figure. The openings, approximately 





¢ Pes pg 
Fig. 1. Cross section of wood (1000) showing an aspirated bordered 
pit connecting two adjacent fiber cavities. 


circular in shape, in the adjacent cell walls are separated by a membrane that 
is a continuation of the middle lamella—the cementing material between the 
fibers. This membrane, particularly with softwoods, has a slight central thick- 
ening called the torus. In some cases the pits are aspirated, that is the torus 
is held against one of the seats of the cell wall opening, thus making the pit 
act like a closed valve (Fig. 1). 

The permeability of a pit will depend upon its condition. If the pit is not 
aspirated, the permeability will be a function of the porosity, the area, and 
the thickness of the membrane. If it is aspirated, the permeability will fur- 
ther depend upon how closely the valve seat fits and whether it is cemented 
in a closed position with resinous material. 

Consideration of the fact that Fig. 1 is magnified to the limit of reasonable 
definition makes obvious the additional fact that an examination of photo- 
micrographs will give little information on permeability. Such a method of 
study is eliminated also by the large number of measurements that would 
have to be made. Each softwood fiber with an average length of about 0.3 cm 
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and a diameter of approximately one hundredth of this value has from 30 to 
300 pits connecting it to adjoining fibers; and there are from 50,000 to 100,000 
of such fibers in a square centimeter of cross section. 

The fine capillary structure of wood, however, can be fairly well char- 
acterized by a combination of three different types of physical measurements. 
In the previously reported researches from this Laboratory on capillary struc- 
ture of wood, the effective continuous capillary cross section was determined 
indirectly by means of electroendosmotic flow measurements.?* This method 
has recently been abandoned because it is involved and requires several ap- 
proximations. The ratio of the effective capillary length to the effective con- 
tinuous capillary cross section is now being determined by measuring the elec- 
trical resistance of strong salt solutions in the wood structure and in bulk.*® 
An average effective capillary radius (fourth root of the average fourth 
power) is determined by combining the effective capillary values with data 
obtained by pressure-permeability measurements of a continuous water phase 
extending through the wood sections, using Poiseuille’s law.?* The maximum 
effective capillary radius is determined by measuring the gas pressure re- 
quired to overcome the effect of the surface tension of water in the capillary 
systems.?* An additional method of study, consisting of measurements of the 
permeability of the wood to colloidal solutions containing particles of ap- 
proximately known size, serves as a check upon the capillary sizes obtained 
by the two previous methods. Since no recent measurements by this method 
have been made, the method will be dismissed here with the statement that 
it has always given results consistent with those of the other methods.?*:® 

The data presented in this paper will be confined to the simplest case, 
coniferous woods and transverse sections in which the flow of electric cur- 
rent, flow of water, and displacement of water are all in the direction of the 
fiber length. The specimens were all cut from fresh, unseasoned wood. 


DETERMINATION OF THE Ratio or THE EFFECTIVE CAPILLARY 
LENGTH TO THE ErrectivE Continuous CAPILLARY Cross 
SECTION BY AN ELEectrIcAL Conpuctivity MEtTHopD 

The ratio of the effective capillary length to the effective continuous capil- 
lary cross section can be calculated from the electrical resistance of the wood 
sections, the voids of which are completely filled with a salt solution, and from 
the specific resistance of the salt solution in bulk, providing the surface con- 
ductivity is made negligible and the conductance of the cell wall is negligible. 
The surface conductivity for very dilute salt solutions and especially for dis- 
tilled water may be many times the bulk conductivity of the solution or of the 


2 Stamm, Colloid Symposium Monograph 6, 83 (1928). 

* Stamm, J. Agr. Research 38, 23 (1929). 

4 Fairbrother and Mastin, J. Chem. Soc. 125, 2319 (1924). 

5 Marshall, J. Soc. Chem. Ind. 48, 373T (1929). 

6 Bechold & Szidon, Kolloid Z., Special No. April 1, 259 (1925); McBain and Kistler, J. 
Phys. Chem. 35, 130 (1931). 
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water,?:45.7.5 but for appreciable salt solution concentrations this surface con- 
ductivity becomes negligible in comparison with the bulk conductivity. When 
the potassium chloride solutions used in this investigation exceeded a con- 
centration of 0.07 mols per 1, the ratio of the specific resistance of the salt 
solution in bulk to its resistance in the wood structure was found to be inde- 
pendent of concentration, thus indicating that the surface effects are negligi- 
ble. Hence, measurements were made using salt solutions exceeding this con- 


centration. 

The specific conductance of dry wood is extremely small, approximately 
10-” mho on an over-all dimension basis. When water is adsorbed by the cell 
wall the conductivity increases so as to give a linear relationship between the 
logarithm of the conductivity and the moisture content, up to fiber satura- 
tion, at which point the specific conductance is 310~-* to 110-5 mho.® 
The water adsorbed in the cell wall at this point shows an increased conduc- 
tivity over the bulk conductivity of water, because of surface conduction. 
Salt solutions, however, act differently. The conductivity of such solutions 





Fig. 2. Diagramatic representation of the capillary path through wood showing the manner 
in which the fiber cavities are connected through the pores of the pit membranes. 


at the fiber-saturation point, although exceeding that for water, is sufficiently 
less than the bulk conductivity to make the conductivity of the cell wall in 
the presence of the salt negligible for the present measurements. This is due 
to the fact that salt solutions are not dispersed in the cell wall substance in 
a continuous manner as with water, but are confined to the grosser capillary 
spaces; this fact will be considered more at length in a later publication. 

The electrical conductivity of a transverse wood section, the capillary 
structure of which is filled with a salt solution, is thus substantially equal to 
the sum of the bulk conductivities of the solution in all of the individual capil- 
lary paths connected in parallel. These capillary paths in turn are made up 
of fiber cavities and pores of pit membranes connected in series (Fig. 2). 
The part of the electrical resistance for which the fiber cavities are responsible 
can be calculated from the fractional void volume or average fractional capil- 
lary cross section of soaked swollen wood, g, which in turn can be obtained 
from the bulk density d of the wood on a wet volume and dry weight basis. 
and the density of wood substance, dp which is equal to 1.52. 

Thus, 


1 : (1) 
[7-2 
do 
7D. R. Briggs, J. Phys. Chem. 32, 641 (1928). 
® McBain, Peaker and King, J. Am. Chem. Soc. 51, 3294 (1929); McBain and Peaker, J. 
Phys. Chem. 34, 1033 (1930). 
® Stamm, Ind. Eng. Chem, Anal. 1, 94 (1929). 
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The value of g obtained in this way includes the voids in the cell wall, which 
as already shown are electrically ineffective, as well as the fiber cavity voids 
and other void structure which, for this purpose, is negligible. The fractional 
void volume or average fractional void cross section of the cell wall is equal 
to the product of the moisture content M per gram of dry wood at the fiber- 
saturation point, and the density of the wood, d. The average fractional cross 
section of the fiber cavities, 


qe =q— Md (2) 
and the electrical resistance of the fiber cavities 
Raf, 
Ry = (3) 





(q — Md)Q 


where R,, is the specific resistance of the salt solution, ZL the thickness of the 
section under investigation, and Q the cross section of the specimen. 

The total resistance R,, of the salt solution in the pores of the pit mem- 
branes will depend upon the fractional cross section of such pores traversed 
in parallel and the length of path, that is, the continuous effective capillary 
cross section gm of the pores expressed as a fraction of the cross section of the 





Fig. 3. Cells used in electrical resistance measurements. 


specimen, and the sum /,, of the thicknesses of the pit membranes traversed 
in series. Then, 


R,>l m 
Qgm 


The experimentally measured resistance R is equal to the sum of the resis- 
tances from Eqs. (3) and (4). 





(4) 


™m 











Rs pL Rsglm 
- : (5) 
(q — Md)Q Ogm 
and 
b. ) Rap 
fa O(p__Rot_) 2 
Qu Rep (gq — Md)Q 


The measurements were made upon cylindrical sections of wood that had 
been turned in a lathe in a soaked and swollen condition to fit tightly into 
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the ends of two hard-rubber cells (Fig. 3). Two sets of cells and electrodes 
for sections of 0.665 cm? and 2.66 cm? cross-sectional area were used. The 
electrodes E are heavy platinum wire disk coils sealed into glass tubes T that 
are filled with mercury for electrical contact. The electrodes were freshly 
platinized each day before using. The side tubes S are for convenience in filling 
with liquid. The measuring apparatus consisted of a student circular slide- 
wire bridge, a 4-dial resistance box (0.1 to 999.9 ohm), a microphone hummer, 
telephone receivers, a condenser, a switch, and dry cells connected according 
to standard conductivity measuring practice. 

The sections were soaked in distilled water for at least two weeks. To 
facilitate the replacement of air by water the soaking was done in a vacuum 
desiccator to which suction was applied intermittently. Potassium chloride 
was then added. Electrical resistance measurements on the wood sections 
showed that diffusion was complete and equilibrium of the salt distribution 
was obtained in less than a week. The specific resistance of the salt solution 
in equilibrium with the wood sections was determined, using the hard-rubber 
cells of Figure 3 clamped together with a ring gasket replacing the wood sec- 
tion. Measurements were made with the electrodes separated by different 
distances. The measurements were made at 25°C in a thermostatic air bath. 
TABLE I. Ratio of the effective capillary length to the effective continuous capillary cross section 


obtained from electrical resistance measurements at 25°C. KCl conc. =0.79N, 
Rs = 9.84 am, M=0. 30. 














Speci imen | I q Q | L R R; Ln /Qrn 
Slash pine 
Heartwood | 0.430 | 0.717 2. 6th | 1.965 | 15.30 12.40 0.783 
| 0.430 0.717 2.66 1.450 | 11.30 9.13 | 0.587 
0.430 0.717 2.66 | .922 | 7.20 5.81 0.375 
0.430 0.717 2.66 .585 4.68 | 3.69 0.267 
0.430 0.717 .665 1.925 | 59.90 | 48.50 0.771 
| 0.430 0.717 .665 | 1.400 | 43.80 35.30 0.575 
Slash pine | 
Sapwood 0.456 0.700 2.66 1.940 | 15.00 | 12.75 0.608 
| 0.456 0.700 2.66 1.365 | 10.65 | 8.97 0.453 
0.456 0.700 2.66 .960 | 7.38 | 6.30 0.292 
yor 0.700 2.66 .640 | 5.05 | 4.20 | 0.229 
0.456 0.700 .665 | 1.890 58.50 | 49.80 | 0.588 
0.456 0.700 .656 | 1.395 43 .00 | 36.70 | 


rhe data for both the sapwood and the heartwood of slash pine are given 
in Table I and plotted in Fig. 4. Within the range of experimental error the 
straight lines pass through the origin as would be expected according to 
theory. The linear relationship indicates that g, the continuous effective 
capillary cross section, varies but slightly for the adjoining sections and that 
lm, the total length of the effective capillaries, varies directly with the thick- 
ness of the sections, thus indicating a rather uniform distribution of the effec- 
tive capillaries. Such a relationship is rather to be expected, considering the 
larze number of fiber cavities and pit membranes traversed by the current. 
It is rather surprising, however, that the difference between the sapwood and 
the heartwood is-so small. This fact can be better understood after considera- 
tion of the other capillary methods. 
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Fig. 4. Relationship between the ratio of the effective capillary length to the effec- 
tive capillary cross section and the thickness of the section. 


An approximation of ‘the actual value of g, can be made by taking the 
value of /m/qm from the graphs for a section whose thickness is equal to that 
of half the average fiber length (0.175 cm) and taking the thickness of a single 
pit membrane, which was found microscopically to range from 2X10- to 
210-4 cm. Dividing Im by Jm/qm gives a value of gm ranging from 3.61074 
to 3.6X10-*. This is in as good agreement with the value of gm obtained by 
electroendosmosis** as could be expected, because of the uncertainty of the 
effective value of /,,. 


DETERMINATION OF THE AVERAGE EFFECTIVE CAPILLARY SIZE BY 
Means or Hyprostatic FLow StupiEs 


In the capillary flow of liquids under hydrostatic pressure the velocity 
is dependent not only upon the total continuous effective capillary cross sec- 
tion but also upon the effective size of the individual ducts. By combining 
the previous data for the ratio of the effective capillary length to the effective 
continuous capillary cross section with data from hydrostatic-flow studies, 
it becomes possible to calculate the actual size of the average effective capil- 
lary openings. The passage of liquids under external pressure through a smooth 
capillary tube of circular section for conditions of stream-line flow (which are 
obtained for the capillary systems considered, a correction for impact turbu- 
lence alone being necessary**) may be quantitatively expressed by Poiseuille’s 
law?:3.19 
m'P ArP 


=— = — 7 
8nl 8nl ) 


+ 





in which V represents the rate of flow, 7 the capillary radius, / the capillary 
length, P the applied pressure, 7 the viscosity of the liquid, and A the cross- 
10 Reynolds, Roy. Soc. (London) Phil. Trans. 174, 935 (1883); Biglow, J. Amer. Chem. Soc. 


29, 1675 (1907); Bartell, J. Phys. Chem. 16, 318 (1912); Bartell and Carpenter, J. Phys. Chem. 
27, 252 (1923); Hitchcock, J. Gen. Physiol. 9, 755 (1926). 
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sectional capillary area, all in c.g.s. units. For a bundle of N capillaries Eq. (7) 

becomes, 

_  NartP  Qqr’P 
8nl 8yl 


where Q is the total cross section of the bundle of capillary tubes and q is 
the fractional open cross section. In this case r is the fourth root of the average 
fourth power of all of the radii rather than an arithmetic average. Since the 
range in size of the pores of the pit membranes is not great as explained in 
the following section of this article, the value will not differ appreciably from 
the arithmetic average. 

Instead of determining the volume of discharge of water from the capil- 
lary system per unit of time under a definite head, with the system held at a 
constant temperature so that 7 is both determinable and constant, a differ- 
ential method, which involves several special features, was used. If the sec- 
tion under test and a calibrated capillary tube are connected in series with 


Bi OqriPi 


J —_—_——_ 9 
1 Sal; (9) 


/ 





(8) 


for a test section where qi, 71, and /,; represent the effective fractional cross 
section, effective radius and effective length, and 


en. (10) 
vl Syl 
for a calibrated standard capillary tube then since 
V; = Ve (11) 
mrs'l, Ps 1/2 
nn (aR) va 
Qleqi Pi 


The question arises as to what constitutes the effective cross section, effec- 
tive radius, and effective length for hydrostatic flow through the broken capil- 
lary system of wood. This question is answered by a simple experiment of 
Ewert" who found that under a head of 80 cm of water 4.0 to 4.2 cc of water 
flowed per hour through a capillary tube having a radius of 0.0075 cm and a 
length of 8 cm; the theoretical flow for these conditions was 4.3 cc per hour. 
The tube was then cut into 16 lengths of 0.5 cm and each length was sealed 
into a tube 0.2 cm in radius and 5 cm long. All of these units were connected 
in series, thus giving a discontinuous capillary analogous to the capillary path 
through a section of wood. Under the same head of water, 3.8 to 3.9 cc passed 
through the tube per hour, showing substantially that only the fine capillaries 
are effective in determining the rate of flow and that the rate is practically 
the same as for one capillary with a length equal to the sum of the lengths of 
the sections, the small variation being accountable for on the basis of increased 
impact turbulence, which is corrected for in this investigation. 

When the frictional resistance in the fiber cavities is negligible in com- 


11 Ewart, Roy. Soc. (London) Phil. Trans. (B) 198, 41 (1905). 
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parison with that in the pores of the pit membranes, as it was in the preceding 
illustration, the effective dimensions q:, 71, and /; are equal to the correspond- 
ing dimensions gm, fm, and 1, of the pores of the pit membranes. Eq. (10) then 


becomes 
Tr2'lmP 2 ass 
_n = <=) , (13) 
QlegmP 1 


Where the resistance to flow through the fiber cavities is not negligible in 
comparison with the resistance through the pit membranes, Eq. (10) assumes 
a more complicated form 




















1 1 1 
——e | Sp . . ae . (14) 
tre'P» Ogyr Pi Ogmrn*P 
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and 











Im i — 

rm = (- \ - (15) 

Ym ( Ql,P L ) 
tr2'P» gsr s” 


where g; and r,; are for the fiber cavity dimensions, and L is the thickness of 
the section. g; is equal to (¢— Md) as in Eq. (2) and 7; is the average radius 
of the fiber cavity estimated from Table 10 of the previous paper.* 

The apparatus used for these measurements is shown in Fig. 5. It consists 
essentially of three tubular glass compartments, I, II, and III, separated by 
means of the wood section M and the standard capillary tube C. A new clamp- 
ing device D is made of brass and is fastened to tubes I and II with de 
Khotinsky cement. Rubber gaskets are used on each of the clamp faces. This 
clamping device not only has the advantage over the previously used type in 
forming a tighter seal, but also in preventing evaporation from the edges of 
the section. Mercury manometers P; and P, indicate the pressure drops be- 
tween compartments I and II; and II and III, respectively. Where the resis- 
tance to flow through the standard glass capillary is very small as compared 
to that through the wood section, thus giving a very small pressure drop 
across P:, the manometer P: is replaced by an inverted water manometer, 
thus increasing the accuracy of reading 12.6 times. A is the intake and B the 
over-flow. The tubes S are merely for convenience in filling the apparatus 
with water and displacing the air. A constant hydrostatic head of water is 
applied at A. When the pressure drops through the standard capillary and 
the wood section came to equilibrium, which required several hours for resist- 
ant heartwood sections but only a few minutes for the sapwood sections 
readings on the manometer tubes were recorded. Reversing the connections 
so that the hydrostatic head is applied at B on the standard capillary side 
rather than on the wood-section side, with A serving as the over-flow, has 
no effect upon the equilibrium pressure-drop ratio. 

The pressure-drop ratio P,/P, is not constant for different applied heads, 
as would be expected for an ideal system free from turbulence. This deviation 
from constancy was shown to be due to impact turbulence at the ends of the 
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capillary tubes.?’* The impact turbulence is a function of the velocity of flow 
and of the length of the capillary (see Couette’s amplification of Poiseuille’s 
equation”) and becomes negligible for infinitesimal applied heads. Hence, it 
can be eliminated by determining P:/P, for a number of different applied 
heads and extrapolating to zero applied head. This correction is explained 
more in detail in an earlier publication.” 


























Fig. 5. Differential pressure-drop hydrostatic-flow apparatus. 


For these measurements it is extremely important to displace all of the 
air from the sections. This was done by intermittently applying suction in a 
vacuum desiccator to the sections immersed in water. Heating the water to 
60°C facilitated the removal of air. This treatment in all cases was continued 
for at least two weeks. 


TABLE II. Effective capillary sizes from hydrostatic flow and from the overcoming of surface tension, 
Q=2.08 cm? (for hydrostatic flow), 77 =0.001 cm. 


























| | | Im/qm | Ps/Ps | | 
Specimen | L | rot/l,z |fromFig.| (extra- | te | p aaa 
| | + | polated) | 
Cm Kg/cm? Cm : 
Slash pine 1.340 (0.53X10-%} 0.54 0.013 2.37 X10-4| 27.0 5.5 X 10-8 
Heartwood 0.875 53 0.35 0.025 2.65 22.0 6.7 
.670 | .53 0.27 0.037 2.83 15.5 9.5 
.670 |3.73 0.27 | 0.0047 | 2.66 —_ on 
| Cm of 
| mercury 
Slash pine 1.980 |32.1K10- 0.62 1.25 2.19X10-4} 12.0 9.0 10-4 
Sapwood 1.365 {32.1 0.43 1.70 | 2.11 11.2 9.7 
0.915 /32.1 0.29 2.40 2.04 10.2 10.6 
.730 132.1 0.23 2.85 | 1.97 9.5 11.4 

















2 Herschel, Colloid Chemistry, Theoretical and Applied 1, 727 (1926) edited by Alexander 
Chemical Catalog Co., New York. 
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The average effective radii of the pit-membrane pores of the heartwood 
and the sapwood of slash pine are given in Table II. With heartwood, the pit 
membranes are alone effective in resisting flow, whereas with sapwood the 
fiber cavities are also partially effective. The differences in the effective pore 
sizes for these two parts of the same wood are rather large. 


DETERMINATION OF THE Maximum Errective CAPILLARY SIZE BY 
OVERCOMING THE Errect oF SURFACE TENSION WITH 
APPLIED GAS PRESSURE 


The determination of the pressure required to overcome the effect of sur- 
face tension in the capillary system permits the calculation of the maximum 
effective capillary diameters, using Jurin’s law. 

Liquids with a density d that wet capillary tubes with an angle of wetting 
6 rise in the tubes to a height dependent upon the radius of the tube ,;, 
the surface tension of the liquid a, and the counter pull of gravity g. Equi- 
librium is reached when the resultant force F is equal to zero. Hence, 


F = rr’hdg — 2xra cos 6 = 0 (16) 
and 
rhdg 


c= 
2 cos 6 





(17) 


For water and other liquids that wet the capillaries completely, # is equal to 
zero and the expression simplifies to, 


rhdg 


¢= “-— (18) 


When @ is known for the liquid used, this equation serves as a means of de- 
termining r by merely measuring /. This method, however, can not be used 
with wood capillaries because their opacity prevents the measurement of h 
and also because the irregularity of the capillary systems makes it practically 
impossible to bring the menisci to the points desired for measurement. Al- 
though h is not directly measureable, the pressure p that will just cause the 
removal of the liquid from a capillary, thus overcoming the effect of surface 
tension, can be readily determined." Since 


p = hdg (19) 
then 
2e 
yu: (20) 
p 


If several capillaries of different bores are combined in a bundle, water will 
be displaced first from the tube of least capillary resistance, which will be 
the tube having the maximum effective bore, the effective bore being the mini- 
mum bore for a capillary of nonuniform cross section. 


‘8 Biglow and Bartell, J. Ame#. Chem. Soc. 31, 1194 (1909); Bartell, J. Phys. Chem. 16, 
318 (1912); Bartell and Carpenter, J. Phys. Chem. 27, 252 (1923). 
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The apparatus used for these measurements is shown in Fig. 6. It consists 
of a T tube made of brass, with a flange on one arm for clamping wood sec- 
tions securely over the open end of the arm. The disk plate that holds the 
wood section in place by means of four screws has in it several small holes 
distributed over the area of the tube to allow the free passage of air, while 
still supporting the section sufficiently to prevent rupture from taking place 
readily at high pressures. Rubber gaskets are used between the flange and 
the wood section and between the plate and the wood section so as to make 
the seal practically air-tight. For wood sections of low resistance the gas pres- 
sure was obtained from the laboratory compressed-air supply line and the 
pressure applied was measured on a mercury manometer. For resistant sec- 
tions a tank of oxygen under high pressure, having a readily controllable 
needle valve, formed the source of gas pressure and the pressure applied was 
measured by means of a calibrated pressure gauge, which was sensitive to 
0.5 Kg. 
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Fig. 6. Apparatus for overcoming the effect of surface tension. 


The measurements were made on water-soaked sections that were clamped 
securely between the flange and the plate, and dipped into a dish of water 
just far enough for the liquid to touch the bottom of the wood section. The 
gas pressure was applied gradually, so as to allow for complete dissipation, 
and the pressure at which the first gas bubbles emerged from under the im- 
mersed plate was observed. It is readily possible to distinguish actual passage 
of gas through the section from the slight gas leak between the non-immersed 
flange and the wood; Jeaks of this nature can not be entirely eliminated at 
higher pressures. The possibility of distinguishing the two kinds of gas flow 
makes this type of flange connection preferable to the cup type used with the 
previously described apparatus. After making a determination the apparatus 
was disconnected from the pressure line and a drop of water was placed in 
the tube to assure the resoaking of the upper part of the section; then after 
standing for a few minutes a check determination was made. 

The right-hand part of Table II gives the results of these measurements, 
made at approximately 25°C, for the same heartwood and sapwood sections 
of slash pine used for the hydrostatic flow measurements. The value for the 
maximum effective pore radius in each case decreases with an increase in 
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thickness of the section. This is to be expected; the most effective path 
through a number of pit membranes in series will approach the average effec- 
tive path as the number of pit membranes in series increases, because of the 
decreasing probability of all the pit membranes in series containing pores 
of maximum size. Measurements made by this last method of study upon 
radial sections, in which the tangential displacement of water by gas is 
through the same structure except that far more pit membranes are traversed 
in series, gave values of r for the heartwood of 4.0 10~-° cm and for the sap- 
wood of 3.710-* cm. The number of pit membranes traversed in series for 
these sections was approximately 40, whereas the number traversed for the 
sections of Table II was from 2 to 5. 


CONCLUSIONS 

Data obtained by the three methods used for studying the capillary struc- 
ture of wood indicate that there is a large difference in permeability between 
the heartwood and the sapwood of slash pine; that of the sapwood is 4,000 
to 5,000 times that of the heartwood. This difference in permeability is not 
due to a difference in the effective continuous capillary cross section but 
rather to the difference in the size of the effective openings. The radii of the 
maximum effective openings range from 3 to 6 times the corresponding aver- 
age values. The values for the size of the effective openings obtained by the 
method of overcoming the effect of surface tension approach the average 
values in the case of thick sections where the number of pit membrnaes tra- 
versed in series is large. The fact that two entirely different methods of meas- 
uring the size of the effective openings give results of the same order of magni- 
tude serves as a check upon the validity of the methods. 


SUMMARY 


1. Three physical methods for studying capillary structure have been de- 
veloped and applied to the study of the fine structure of wood. 

2. The ratio of the effective capillary length to the effective continuous 
capillary cross section is obtained from the electrical resistance of salt solu- 
tions filling the wood structure and the specific resistance of the solutions in 
bulk. 

3. The average effective capillary radius is obtained combining data from 
the preceding measurements with data obtained from hydrostatic-flow studies 
of a continuous liquid phase extending through the wood sections, using 
Poiseuille’s law. A differential pressure-drop apparatus was devised for these 
measurements, in which the flow of liquid through a standard capillary tube 
and the test specimen connected in series is determined by measuring the 
pressure drop across each. 

4. The maximum effective capillary radius is determined by measuring 
the gas pressure required to overcome the effect of the surface tension of a 
liquid in the capillary systems. 

5. Data obtained by these methods show the large difference between the 
fine capillary structure of heartwood and of sapwood. 
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ABSTRACT 


A study of the swelling and gelling properties of bentonite dispersions has been 
made in connection with the increasing use of this mineral for drilling oil wells. It is 
shown that the emulsoid type of colloids present responds to changes in the pH of the 
dispersion medium with respect to rate of settling, viscosity and swelling, thus behav- 
ing as a typical hydrophilic colloid. The curve for percent of solid matter settled in a 
definite time as ordinate against the pH shows a maximum at the isoelectric point, 
with minima on each side of this point, corresponding to the points of maximum ad- 
sorption. As the pH is further increased or decreased flocculation begins with increased 
rate of settling. The viscosity curve is closely related to the rate of settling curve, 
showing maxima and minima corresponding to the electro-viscous effect. Maximum 
swelling occurs at approximately the point of maximum viscosity. Methods are given 
for controlling the swelling of bentonites. Relaxation curves are given for bentonite 
suspensions and utilized to explain the suspension qualities of these sols for heavy min- 
eral dispersions, such as barytes or amorphous silica. A few applications of bentonite 
sols to the drilling of oil wells are given. 


1. INTRODUCTION 


HE examination of the properties of bentonite suspensions is of interest 

on account of its swelling and gelling properties, together with its multi- 
tude of industrial uses.' In this connection our particular interest is the use of 
the mineral in rotary drilling fluids for sinking oil wells. By virtue of their 
colloidal properties natural and prepared bentonites are having an increased 
application in the technique of oil-field drilling, both as protective colloids 
for the stabilization of suspensions and alone in water as a self-contained 
drilling fluid. Any drilling fluid may be considered as a three-component sys- 
tem: water, gel-forming colloids, and non-gel-forming colloids. Of these two 
types of colloids the gel-forming type is by far the more important, for this 
group is responsible for the necessary characteristics of a drilling fluid, i.e. 
viscosity, suspension, body or consistency, lubrication, sealing and gel proper- 
ties in general. Colloids of the non-gelling type, together with larger particles, 
contribute partially to the viscosity but more particularly to the specific grav- 
ity of the drilling fluid. The understanding of the proper use of bentonite in 
drilling fluids demands a knowledge of the properties of the gel-forming col- 
loids in particular. We have, therefore, examined natural and treated bento- 


1 Davis and Vacher, Bur. Mines Tech. Paper 438; Alexander, Colloid Symposium Mono- 
graph 2, 99 (1925). 
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nite for the purpose of gaining a clearer comprehension of the physical and 
chemical properties of this mineral with respect to its emulsoid content and 
how this type of colloid may be chemically and physically controlled to have 
maximum efficiency for the practical use under consideration. 


Il. PRoPERTIES OF BENTONITE 


Certain bentonites of high colloidal content have astonishing powers of 
imbibition. A 5 percent dispersion in water will swell and form a stiff gel, oc- 
cupying a volume twenty times the original volume of the solid mineral. Small 
percentages of bentonite added to clay dispersions that have lost their quality 
of low settling rate, due to loss of colloidal material, will restore these quali- 
ties, and the viscosity of the suspension is considerably increased. This finds 
practical application in restoring spent clay drilling fluids that have lost their 
colloids in porous formations. A clay dispersion requires about ten times as 
much weight of solid material in suspension as a bentonite dispersion to have 
the same viscosity as the latter. 

In drilling through formations of high gas pressure weighted drilling fluids 
are necessary to furnish the required hydrostatic head to balance this pres- 
sure. Heavy minerals such as barytes are used for this purpose, but the sus- 
pension qualities of such materials are poor. A small percentage of bentonite 
added to such suspensions will entirely arrest the tendency to settle by virtue 
of the protective action of the bentonite colloids. 


Ill. Origin AND OccURRENCE OF BENTONITE 


The principal commercial deposits of bentonite with high content of col- 
loids are in South Dakota, Wyoming and California. It appears that the 
mineral is a decomposition product of volcanic ash.? There are several reasons 
comprising this evidence, such as the absence of quartz, which is present in 
all ordinary clay sediments, the presence of feldspar, and the extension of 
individual beds over great distances without much variation in thickness. 
The clayey matter in normal sediments has been so acted on by electrolytes 
and other active substances that it has lost its original porosity and tendency 
to imbibe water with attendent swelling, which bentonite exhibits to such a 
marked degree. The chemical composition of bentonite is, moreover, exactly 
what one would expect from the partial alteration of an intermediate igneous 
rock, accompanied by the removal of alkalies and addition of water. 

Bentonites are probably hydrated silicates of aluminum, according to 
Davis and Vacher,* and may be represented by the formula Al,O3-4SiO,- 
xX H:O. Some bentonites appear to contain about 50 percent of particles of 
colloidal size (5-200my). According to Wherry‘ the colloidal properties of 
bentonite are due to a “felted” texture, the particles crystalline, of colloidal 
thickness, but of appreciable size in the other two dimensions. The adsorp- 
tive capacity for water is supposedly due to the wetting of these flat surfaces 


2 Hewett, J. Wash. Acad. Sci. 7, 196 (1917); Wherry, ibid. 576, 
3 Reference 1. 
* Wherry, Am. Mineral. 10, 120 (1925). 
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by adsorbed water. The properties of gel formation by bentonite suspensions 
is probably due to the interlocking of crystals when the suspension is at rest. 


IV. ExperRIMENTAL WorK 
1. Measurements of rate of settling and viscosity as functions of the p //. 


The colloidal properties of a Wyoming bentonite were examined as a func- 
tion of the acidity or alkalinity of its water dispersion. The pH was measured 
by means of the hydrogen electrode and calomel cell with alkaline dispersions 
and a quinhydrone electrode and calomel cell with acid dispersions. By meas- 
uring the pH of the dispersion and the aqueous phase separated from it by 
centrifuging, it was found that they were identical, thus proving that equi- 
librium was attained in this adsorption process. 

The stability of the suspension may be measured by the rate of settling 
of solids when centrifuged at high speed. Suspensions of bentonite containing 
5 percent by weight of solids were centrifuged for 5 hours, after which the per- 
cent settled and the pH were measured. Fig. 1 shows the percentage settled 
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Fig. 1. Percent settled in 5 hours vs. pH. 
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in 5 hours as a function of the pH. A colloid is most stable (deflocculated) 
at the point of minimum rate of settling. The curve shows minimum rates of 
settling both in acid and in alkaline solutions. The colloid is in its defloc- 
culated condition between pH =5.8 and pH = 11.4. At values of pH beyond 
these limits floccuation is attended by increased rates of settling. 

The fact that a bentonite swells in water indicates the presence of a cer- 
tain percentage of colloids of the emulsoid type, stabilized by adsorption of 
salts contained in the mineral. A certain percentage of the colloids are of the 
suspensoid type, and these are in excess. Fig. 1 refers to the behavior of the 
emulsoid content of the suspension. Rate of settling in the deflocculated re- 
gion is at a maximum at the isoelectric point, of course, for stability depends 
upon the charge on the particles and the viscosity. As the pH increases from 
the isoelectric point, the hydroxy] ion is increasingly adsorbed, with the rate 
of settling reaching a minimum value at the limit of deflocculation on the 
alkaline side, which is also the point of maximum adsorption of the hydroxyl 
ion. In more strongly alkaline solutions the cation begins a discharging effect 
on the negatively-charged emulsoid and flocculation begins, with increase in 
the rate of settling. Similar relations hold on the acid side of the isoelectric 
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point. In such a colloidal suspension rates of settling are little dependent upon 
Stokes’ law governing the rate of fall of particles in a liquid, for the lower 
limit of particle size to which this law applies is dependent not only upon 
viscosities of the fluids but more upon the state of dispersion of the particles. 
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Fig. 2. Viscosity vs. pH. 


It is evident from a comparison with Fig. 2 that the viscosities measured 
dynamically are not sufficient to account for the large changes in rates of 
settling as the pH is varied, although there is a close relation between the 
two curves. 

In Fig. 2 it is seen that a minimum viscosity occurs at approximately 
the pH of maximum stabilization as shown in Fig. 1 at pH =11.4. This is to 
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Fig. 3. pH of bentonite suspensions vs. alkali added. 


be expected, for the electro-viscous effect is necessarily a minimum at this 
point. On the addition of acid or alkali, attended by adsorption and increase 
of positive or negative charge, respectively, an increased viscosity is ex- 
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pected in accordance with Smoluchowski’s®’ modified Einstein equation for 
the effect of electric charges of particles on the viscosity of the fluid. This 
equation for the electro-viscous effect has been experimentally verified by 
Kruyt.® As the maximum charge is dissipated on further increase in the pH, 
the viscosity must again decrease, reaching a minimum and rising again very 
rapidly when the colloid is flocculated with gel formation, as shown by the 
curve on the alkaline side. 

The action of alkali upon the bentonite is evidently adsorption. The meas- 
ured pH is not a measure of the alkali added, for alkali disappears from solu- 
tion. This adsorption of alkali by bentonite is shown in Fig. 3, where pi is 
plotted against added alkali in terms of amount of caustic soda in mols per 
liter, giving a typical adsorption isotherm. 

The viscosity curves of the suspensions indicate that no absolute measure 
of viscosity can be made of a gelling colloidal suspension. Fig. 4 shows the 
measured viscosities of two bentonite suspensions with varying speeds of the 
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Fig. 4. Viscosities of bentonite suspensions vs. revolutions per second of viscometer. 


Stormer viscosimeter; the viscosity is represented as a function of the revolu- 
tions per sec. of the cylinder. The viscosity of each suspension is noted to in- 
crease sharply with decreased speed of the rotating cylinder, i.e. viscosity de- 
creases with increased shear, a general behavior of hydrophilic sols.” This 
appears to indicate that physical forces in the gel oppose the motion of the 
particles in other than viscous flow. The viscosities recorded in Fig. 2 were 
measured at a constant speed of the viscosimeter drum of 10 revolutions pei 
second as a standard. 


2. Measurements of swelling. 


The swelling of bentonite was measured in alkaline solutions to find the 
pH at which maximum swelling would occur. One gram samples of bentonite 
placed in uniform tubes were treated with solutions varying in concentration 


5 Smoluchowski, Kolloid-Z. 18, 194 (1916). 

6 Kruyt, ibid. 31, 338 (1922). 

7 Hess, Kolloid-Z., 27, 154 (1920); Hatschek, ibid. 13, 88 (1913); 7, 301 (1910); 8, 34 
(1911); 11, 280, 284 (1912); 12, 238 (1913). 
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of caustic soda. These were allowed to stand for 48 hours with occasional 
agitation to allow contact between various portions of the liquid and solid. 
The measured volume of the bentonite at the end of this period is plotted 
against the pH of the aqueous phase in Fig. 5. Maximum swelling occurred 
at pH =12.7, approximately the point of maximum viscosity as shown in 
Fig. 2. Flocculation is first attended by swelling and then the swelling in- 
creases with increased alkalinity. A similar maximum is found in the acid 
range, but difficulty was experienced in obtaining the equilibrium pH. Meas- 
urements were not continued in the acid range, for we are interested only in 
the alkaline range insofar as drilling fluids are concerned. 

While it is perfectly possible to control the viscosities and gelling of ben- 
tonite gels by the addition of alkali or lime to the suspension, it is desirable 
from a practical standpoint to treat the solid material in such a way that 
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Fig. 5. Volumes occupied by one gram of bentonite vs. pH. 


rapid gelling takes place on dispersion in water with a minimum amount of 
material and without further adjustment of the pH of the suspension. Most 
bentonites will form a rigid gel in a 10 percent dispersion in water without 
further adjustment of the pH, but a 4-5 percent dispersion will give the same 
results when properly treated. It is interesting to note that we were unable to 
obtain a gel by adding solid caustic soda to the dry bentonite with subsequent 
dispersion of the mixture, although the pH of the resulting dispersion was 
about 12.7. However, addition of 0.5—1 percent magnesium carbonate to the 
dry bentonite will result in a rigid gel in a 5 percent dispersion. No doubt 
the unique properties of magnesium silicates play a role in this case, as well 
as the proper adjustment of the pH. 


V. Errect or SALTS ON THE SWELLING OF BENTONITE 


The effect of salts on bentonite gels has been studied by Davis and 
Vacher.® Such salts as sodium and potassium chlorides show a marked effect 


® Davis and Bacher, reference 1. 
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in prevention of the swelling of bentonites. This effect is particularly desirable 
in the prevention of trouble due to swelling in drilling through bentonite 
formations and materials resembling them. Since bentonites from which the 
salts have been removed by dialysis exhibit neither swelling nor osmotic pres- 
sure, it appears that swelling and water absorption are governed by the solu- 
ble salts associated with the solid material. The use of salts in the external 
phase will naturally decrease the activity of water to be adsorbed on the 
solid surfaces of the mineral crystals, with attendant swelling. 


VI. Use or BENTONITE IN WEIGHTED DRILLING FLUIDS 


One of the most useful properties of bentonite in connection with drilling 
fluids is the stabilization of so-called weighted fluids, suspensions of materials 
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Fig. 6. Amorphous silica suspensions. Density = 1.44. 


such as barite, hematite, amorphous silica, etc. These fluids are used for drill- 
ing through high pressure gas formations, in order to furnish sufficient hydro- 
static pressure to counterbalance the formation pressure. The great difficulty 
with these suspensions is their tendency to settle in the drill hole, making 
further progress in drilling impossible in may cases. The addition of relatively 
small quantities of bentonite to these suspensions has a remarkable stabilizing 
influence, as shown in Fig. 6 for an amorphous silica suspension of density 
1.44; 3 to 4 percent of bentonite (of total water in the suspension) decreases the 
rate of settling of the suspension from approximately 40 percent in 24 hours 
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to practically zero in the same time. The viscosity curve of the suspension 
rises rapidly on the addition of bentonite, as shown in Fig. 6. 

The explanation of this large decrease in rate of settling is furnished by 
the data plotted in Fig. 7, where the deflection in arbitrary units of the sus- 
pended cylinder of a MacMichael viscosimeter is plotted against the time 
in minutes, i.e. the relaxation curve.® A residual deformation is maintained 
permanently due to gel formation of the added bentonite. As seen from the 
curves of Fig. 7 an amorphous suspension of density 1.44 has zero residual 
deformation, while the addition of sufficient bentonite to stabilize the suspen- 
sion results in a definite value for residual deformation. A heavier suspension 
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Fig. 7. Residual torque of amorphous silica suspensions. 


(density 1.64) prepared from amorphous silica alone to have the same viscos- 
ity as the lighter suspension plus bentonite had a very small residual deforma- 
tion value, as shown in Fig. 7 

Gel formation, with attendent arrest of settling of suspensions, by addi- 
tion of bentonite,!° or chemicals in the case of clay suspensions, is finding in- 
creased application in the art of drilling oil wells. 


* See Schwedoff, J. Physique (2) 8, 341 (1899); Hatschek and Jane, Kolloid-Z. 39, 300 
(1926). 

10 Several varieties of treated bentonite are sold for this purpose under various trade- 
names, such as “Aquagel.” 
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THE EDITOR’S COLUMN 


F THE ability to get others to do his work 
] is a mark of anything, your Editor can 
qualify. This month Professors Buchta and 
Hill of the University of Minnesota comment 
in these columns on various topics which have 
been the subject of articles in recent issues of 
the Physical Review. 


The hot summer months are not altogether 
happy ones for initiating anything. I am sorry 
that it will be necessary again to defer the ap- 
pearance of the section “Engineering News 
and notes.” 


em 2 





X-ray Wave-lengths, Crystal Structures and Fundamental Constants 


N ADDITION to discovering new phe- 
] nomena and facts, physicists are attempt- 
ing to describe more accurately those previ- 
ously observed in nature. These descriptions 
are always eventually given in numerical 
terms and the refinements in description in- 
volve increased precision in measurements. In 
the past such increased precision has some- 
times led to the discovery of new phenomena 
and also allowed more critical tests to be made 
of theories relating different phenomena. The 
measurement of x-ray wave-lengths is one 
which can now be made by several methods 
and with high precision. The disagreement in 
the results obtained by using the different 
methods requires that the values of some of 
the fundamental constants be changed from 
their previously accepted values or that some 
new characteristics of crystal structures be 
taken into consideration. 

Some few years ago it was found possible to 
secure x-ray spectra from ruled reflection 
gratings if the glancing angle was very small. 
In the experimental arrangement the beam of 
x-rays strikes the surface at an angle usually 
smaller than one degree. At these small angles 
the K lines from the heavier elements Cr, Cu, 
Mo, etc., undergo total internal reflection. A 
narrow beam of x-rays upon striking the ruled 
grating then yields two or more beams, the 
reflected beam and one or more diffracted 
beams. The angles to be measured are very 


small but by recording the reflected and dif- 
fracted beams on two photographic plates 
spaced a known distance apart and making 
arrangements to record the direct beam on the 
same plates, the measurement of the angles is 
reduced to linear measurements of the spacing 
of lines on the plates and the distance bet ween 
the plates. 

In a paper published in the May 15 issue of 
the PuysicaL Review by J. A. Bearden, the 
results of a long and carefully performed series 
of measurements of the wave-lengths of the 
Cu and Cr K x-rays are given. The method 
used was that described above. Possible errors 
in the grating, those due to divergence of the 
beam, the shrinking of the emulsion on the 
plate, etc. were either eliminated or their 
value determined. The results obtained are 
given in the third column of the following 
table of wave-lengths. 

The wave-lengths as determined by Mr. 
Bearden depend entirely on length measure- 
ments described above and the constant of the 
ruled grating, also a linear measurement. The 
wave-lengths are calculated from the usual 
grating formula 

ny = d(sin zt — sin r) 
or its modification for small angles. In these 
equations n is the order of the spectral line, 
d distance between rulings on the grating and 
z and r the angles of incidence and diffraction 
respectively. There seems to be no reason for 





TABLE I 
Spectral | Crystal | Grating | Limiting | Grating A 
Line nN Error | — Crystal A 
Cu Ka 1.389144 | 1.39225A +0.00014 | +0.224% 
Cu Kae 1.53838 1.54172 +0.00015 +0.217 
Cr Kp | 2.08017 2.08478 | +0.00021 | +0.222 
Cr Kg | 2.28590 2.29097 +0.00023 +0.222 
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questioning the methods of the measurements 
or the use of the measurements in calculating 
the wave-lengths. The wave-lengths so ob- 
tained, however, differ from the previously 
accepted values obtained by using crystals by 
an amount very much larger than the known 
probable errors of either group of measure- 
ments. The differences in percent are given in 
Table I above. The wave-lengths obtained by 
using crystals are consistently shorter than 
those obtained by the use of gratings. 

These differences are in themselves inter- 
esting but are particularly so because of the 
entangling alliances of the crystal wave- 
lengths. The crystal wave-length is deter- 
mined by the Bragg law of diffraction, 


i—»z 
sin? @ 


ny = 2dsiné| 1 — (1) 
in which yw is the index of refraction of the 
x-rays in the crystal, @ the angle of diffrac- 
tion, and d the grating space of the crystal. 
The grating space of the crystal is calculated 
in terms of Avogadro’s number and the crys- 
tal form. Thus for a simple cubic crystal we 
have 

d = (M/Np)' (2) 
in which N is the number of molecules per 
mole, M the molecular weight and p the den- 
sity of the crystal. Using the value of \ from 
ruled grating measurements and putting it in 
equation (1) above we may obtain a value for 
d which when substituted in equation (2) 
gives a value for N, Avogadro’s number. 
When this is done, using constants from 
Birge’s tables, one finds 

N = 6.019 X 10> mol. per mole 
as compared to 
6.064 X 10%, the accepted value. 


Using this value of N and the Faraday con- 
stant F, we find the charge on the electron 
eas 

e = F/N = 4.806 X 10- e.s.u. 
The value given by Millikan is 4.770 x 10-*° 
e.S.U. 

From this value of e we may obtain a corre- 
sponding value of h, Planck’s constant, by 
using the quantum relation 

i-s 


1-—— 
sin? 6 


Ve = hv = hc/X = he/2d sin 0 


in which V is the potential applied to the 
x-ray tube, c the velocity of light and \ the 
minimum wave-length of the continuous x-ray 
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spectrum. Using the results given by Duane, 
Palmer and Yeh for the applied potential V 
and the angle of diffraction we may find the 
value of fA. If in this calculation we use the 
ruled grating wave-length but take e=4.770 
107° e.s.u. we obtain 
h = 6.573 X 1Q°*% erg sec 

as compared to 6.547X10-*’ as given by 
Birge. If we use the ruled grating wave-length 
and the value of e as 4.806107! e.s.u. we 
obtain h=6.623 X 10-?’ erg sec. 

In the first value of h given above we as- 
sumed the accepted value of e rather than the 
value given by grating wave-length measure- 
ments. But even this value of / differs from 
the accepted value by an amount much larger 
than the probable experimental error. 

The values of N, e and the last value of h as 
calculated from the ruled grating wave-lengths 
assume a perfect crystal lattice. It seems at 
present that the only possible way of ac- 
counting for the discrepencies is to assume 
some imperfections or variations in the crystal 
lattice. The constants N, h and e determine 
modern physical calculation through so many 
interlocking directorates that it is difficult 
to accept the high values given by these x-ray 
measurements. Recently Zwicky and others 
have considered variations in the crystal lat- 
tice of crystals to account for their various 
physical properties. In particular, it is as- 
sumed that there is a concentration of atoms 
at rather regularly spaced intervals, giving a 
mosaic structure to the crystal and thus mak- 
ing the average density higher than the den- 
sity of the parts of the crystal effective in 
diffracting x-rays. Such a variation would 
yield crystal wave-lengths that were too short, 
when the calculations assumed a perfect crys- 
tal. The presence of fine cracks in the crystal, 
which may be present in some cases, would 
offer no explanation as the correction neces- 
sary has the wrong sign. 

In conclusion then it seems proper to accept 
Mr. Bearden’s values for the x-ray wave- 
lengths as correct and to assume that the dif- 
ferences of his measurements from the crystal 
wave-lengths are due to crystal imperfections. 
In fact the results may be considered as evi- 
dence for the mosaic structure of crystals. If 
the grating wave-lengths here are correct, the 
true grating space of a calcite crystal is 

d = 3.0359 Anstrom units. 


As this value is 0.22 to 0.25 percent higher 
than values previously used, all wave-lengths 
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based on calcite crystal measurements would 
therefore need to be increased by that frac- 
tion. 


J. A. Bearden, “Absolute Wave-Lengths of 
the Copper and Chromium K-Series,” Phys. 
Rev. 37, 1210 (1931). 

R. T. Birge, “Probable Values of the Gen- 
eral Physical Constants,” Phys. Rev. Supp. 1, 
1 (1929). 
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stants,” Phys. Rev. 35, 1231 (1930). 
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Hyperfine Structure and the Constitution of Nuclei 


structure of the nuclei of the various 
atoms is one of the most weighty questions 
facing physicists today. We now believe that 
there are but two primary elements out of 
which all matter is constructed; i.e. the elec- 
tron and the proton, the latter being the nu- 
cleus of the simplest of the atoms, hydrogen. 
The nuclei of all of the other atoms are sup- 
posed to be made up of electrons and protons. 
From the well-known relations between the 
atomic weights and the atomic numbers of the 
elements as expressed in Mendeleeff’s periodic 
table, one can estimate the number of protons 
and the number of electrons in each of the 
atomic nuclei, but we cannot as yet tell just 
what the dynamical constitution of these sys- 
tems may be. 

Within the last five or six years, however, 
advances have been made in both theory and 
experiment which give promise of ultimate 
success. One of the most interesting of these 
advances has been the interpretation of the 
hyperfine structure of the spectral lines due to 
the atoms. It is well known that under rela- 
tively low dispersion it is quite usually the 
case that the spectral lines of an atom are 
grouped together in sets of a few lines, as for 
example the ordinary yellow doublet of Na at 
5890 and 5896A respectively. This grouping 
of the lines is called a “fine structure” or “mul- 
tiplet structure,” and is quite well accounted 
for by the new quantum mechanics of the 
atom in terms of the structure of the electron 
shells alone. So far as the explanation of this 
feature of the spectrum is concerned, the 
nucleus might just as well be a single point 
charge with charge Ze. If, however, the spec- 
trum is examined with instruments of very 


fie unravelling of the problem of the 


high resolving power it is found that in many, 
although by no means in all, cases the lines of 
the fine structure pattern which before seemed 
single, actually have a complex structure, be- 
ing themselves composed of a series of very 
closely spaced lines. This finer arrangement 
of lines is referred to as a “hyperfine struc- 
ture,” and it is in the explanation of the origin 
of this structure that new light is shed on the 
properties of the nuclei. 

It was suggested in 1927 by W. Pauli, Jr. in 
Germany, that this hyperfine structure might 
be attributed to the presence of an intrinsic 
vector angular momentum and a correspond- 
ing magnetic moment associated with the nu- 
cleus.! This means that the nucleus can no 
longer be treated as a point charge, but ac- 
count must be taken of its complex structure. 
As the electrons fly about in their peregrina- 
tions within the atom they are disturbed 
slightly by the presence of this tiny magnetic 
moment. According to the theory the princi- 
pal result of this interaction is that the energy 
levels which the atom would have if the nu- 
cleus were simply a point charge, are actually 
split up into a number of separate components 
which, however, lie very close together, since 
the electrons are disturbed but very little by 
the magnetic moment of the nucleus. As the 
lines of the spectrum arise from transitions 
between the energy levels of the atom, this 
structure of the energy levels is mirrored in 
the structure of the spectrum, giving rise to 
the hyperfine structure which is observed. 

1 Cf. R. de L. Kronig and S. Frisch, Phys. 
Zeits. 32, 457 (1931), for a summary of the 
present status of the problem of nuclear mo- 
ments for a more general standpoint than was 
possible here. 
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The way in which the lines of the hyperfine 
structure patterns are grouped; i.e. the num- 
ber and separations of the components of 
which each line of the ordinary multiplets are 
composed, depends according to the theory 
upon the kind of motion which the electrons 
perform within the atom. For electron orbits 
of the so-called penetrating s type in which 
the electron is quite near the nucleus for a 
good share of the time, the disturbing effects 
of the nuclear moment are largest and the 
theory is in best agreement with the experi- 
mental results. It was just these cases, of 
course, which found and 
studied first. 

Within the last few months cases in which 
the valence electrons are not in s orbits have 
been examined, and here difficulties have been 
found in the theory. This is emphasized par- 
ticularly in recent papers of Goudsmit? and 
Breit.’ The former examines particularly the 
behavior of the separations of the hyperfine 
structure components for various kinds of 
coupling among the valence electrons, which 
are assumed not to be in s states. It is found 
that the theoretical formulae do not fit the 
experimental data very well particularly for 
Bi and TI, and the conclusion to which Gouds- 
mit comes is that the expression which is used 
for the calculation of the interaction of the 
electron spin with the nuclear spin is some- 
how incorrect and so is to blame for the dis- 
crepancies. 

Breit re-investigates the whole problem 
anew for the particular case of a single val- 
ence electron in a p state; an electron con- 
figuration which gives rise to an ordinary 
doublet term of the type p; and pi, each term 
of which is in turn to show a hyperfine struc- 
ture. Whereas Goudsmit used a non-relativ- 
istic method, Breit bases his calculation on the 
relativistic wave equation of Dirac, finding it 
necessary to use some approximations in order 
to get numerical results since the exact equa- 
tion cannot be readily solved for this problem. 
The final result is a value for the ratio of the 
separation of components of the hyperfine 
structure in the p; term to the separation in 
the pi, term. On comparison with the meas- 
urements of Wulff on TI I this ratio proves to 
be too small by a factor 2, while the non- 
relativistic theory gives a ratio too small by a 


2S. Goudsmit, Phys. Rev. 37, 663 (1931). 
3G. Breit, Phys. Rev. 38, 463 (1931). 


were normally 
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factor 3.4. Thus the new computation im- 
proves the result, but shows that the theory 
is still in need of improvement. 

The work on hyperfine structure gains much 
of its importance from its relation to the 
properties of nuclei. Quite apart from the 
question of the dynamics of the interaction of 
the nucleus and the outer electrons in the 
atom, which is the subject of the papers dis- 
cussed above, the measurement of hyperfine 
patterns is important because it leads to a de- 
termination of the numerical magnitudes of 
the angular momenta of the various nuclei. 
This is of prime importance in the study of 
isotopes, radioactive atoms, and in numerous 
fields. It is fairly certain now that both the 
proton (hydrogen atom nucleus) and the elec- 
tron have themselves intrinsic angular mo- 
menta of amount 3, measured in Bohr units 
h/2m the so-called “spins.” According to the 
theory for the vector composition of the angu- 
lar momenta of electrons and protons in the 
building up of nuclei, it ought to be the case 
that the resultant “spin” angular momentum 
of a nucleus composed of P protons and E 
electrons should be an integer (0, 1, 2, - - -) 
if (P+£) is even, and should be a half-integer 
3, 13, °° +) if (P+£) is odd. All angular mo- 
menta are measured in units h/27, of course. 
Here too there seems to be a lack of agreement 
between present theory and experiment, since 
for the nitrogen nucleus P+E(=14+7 =21) 
is odd, whereas experiment indicates that the 
N nucleus has an integral valued resultant 
angular momentum, probably 1. This par- 
ticular fact was determined not from hyper- 
fine structure measurements, but from Raman 
effect‘ and band spectra. It has been sug- 
gested by Heitler and Herzberg that maybe 
the electron loses its spin when in the nucleus, 
and according to Langer this may come out as 
a result of the new interpretation which has 
been put on quantum mechanics by Schré- 
dinger.® 


E. L. 


‘F. Rasetti, Proc. Nat. Acad. Sci. 
(1929). 

5 R. S. Mulliken, Trans. Faraday Soc. 25, 
634 (1929). 

6 W. Heitler and G. Herzberg, Naturwiss. 
17, 673 (1929). R. M. Langer, Phys. Rev. 
Aug. 15, 1931. E. Schridinger, Sitz. Ber. d. 
Preuss. Akad. 3, 63 (1931). 
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Theories of the Electron 


VER since the latter part of the nine- 
EK, teenth century when the electron the- 
ory began to play a prominent réle in physics, 
mathematical physicists have been trying to 
find a satisfactory model for the electron it- 
self. The calculation of the mass of the elec- 
tron from the properties of its electromagnetic 
field was a high point in the classical theory, 
but was soon quite over-shadowed by Ein- 
stein’s proposal of the relativistic expression 
for the variation of inertial mass with velocity. 
For most purposes the electron has been 
treated simply as a charged mass point. 

When Uhlenbeck and Goudsmit? discovered 
in 1925 that it is very convenient in the in- 
terpretation of the spectra of atoms to as- 
sume that the electron has a mechanical 
angular momentum (“spin”) of amount 3 
(units h/27), the whole question of the con- 
stitution of the electron was re-opened. Since 
that time the conception of electron spin has 
proved its usefulness in many branches of 
atomic physics. Frenkel and Thomas,* using 
the idea of an electron actually in rotation 
about an axis of symmetry, did succeed in cal- 
culating the correct value for the ratio of the 
mechanical angular momentum to the mag- 
netic moment of the electron. Neither of these 
investigators used the ideas of the new quan- 
tum mechanics since these were just being de- 
veloped at that time. Due to the now historic 
work of Dirac,® the theory of the spin has 
actually been incorporated as an integral part 
of quantum mechanics. Dirac’s wave equation 
for the electron has attracted great attention 
since it satisfies the requirements of the special 
theory of relativity and also predicts the ex- 
istence of the spin. Spin and relativity have 
thus become associated in a very intimate 
way. 

When this equation was studied in detail it 
was discovered that it emhodies some predic- 
tions of a very startling sort. Chief among 
these is the prediction that when exposed to 
the influence of electromagnetic radiation, the 

7G. Uhlenbeck and S. Goudsmit, Nature 
117, 264 (1926). 

8 J. Frenkel, Zeits. f. Physik 37, 243 (1926). 
L. Thomas, Phil. Mag. 3, 1 (1927); Nature 
117, 514 (1926). 

* P. A. M. Dirac, Proc. Roy. Soc. A117, 610 
(1928). 


electron ought to act just as though it were 
capable of changing the sign of its charge, and 
so of becoming a particle of charge +e, a 
“positive” electron.!® This was very startling, 
but Dirac! again came to the rescue with the 
idea that perhaps the universe is full of elec- 
trons of this new type, and that when one of 
them changes sign and so becomes an ordinary 
electron, the “hole” which is left in the infinite 
density of “positive” electrons is to be inter- 
preted as being a kind of proton. In many 
ways this is an intriguing theory, but it has 
the obvious difficulty that the mass associated 
with such a “hole” is equal to the mass of the 
electron and not to the mass actually observed 
for the proton. 

Quite recently Schriédinger® has made a 
suggestion which promises to do away with 
this difficulty. This is discussed in a recent 
paper of Langer." The suggestion is of a 
mathematical rather than of a physical na- 
ture, having to do with the restriction of the 
kind of operator allowable in the quantum 
mechanics in such a way that the unpleasant 
transitions of the type discussed above have 
zero probability of occurrence.” As an alter- 
native suggestion, Langer proposes to modify 
the expression for the mutual interaction of 
two particles in such a way that the same re- 
sult is obtained. These theories have not been 
studied in sufficient detail as yet to give a real 
indication of their generality, but they do 
seem to offer possible avenues for escape from 
the difficulties of the earlier work, and are 
indicative of new methods of approach to the 
problem. 


E. L. Beas 


10 Cf. e.g. O. Klein, Zeits. f. Physik 53, 157 
(1929). J. R. Oppenheimer, Phys. Rev. 35, 
461 (1930). 

1 P,. A. M. Dirac, Proc. Roy. Soc. A126, 
360 (1930). 

2 Schrédinger’s paper is not available at 
present to the writer, and is known to him 
only through the discussion of Langer. 

13. It is not clear from Langer’s paper what 
the connection is between Schrédinger’s even 
and odd functions and the symmetry proper- 
ties of the wave equation. It is implied that 
the separation of operators and functions into 
two classes is of a purely arbitrary nature. 

















